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ABSTRACT 


An  investigation  of  the  physical  properties  of  the  superconductive 
lead-indium  alloy  system  has  been  carried  out  as  a  function  of  hydro¬ 
static  pressure,  the  results  of  which  are  presented  in  this  thesis. 

Samples  were  prepared  in  the  form  of  tunnel  junctions  using 
standard  vacuum  deposition  techniques  from  extremely  high  purity 
source  elements.  Pressure  was  applied  to  the  samples  by  suspending 
them  in  a  bath  of  solid  helium  and  their  conductance  characteristics, 
through  which  most  of  the  information  was  gained,  were  recorded  using 
phase  sensitive  harmonic  detection  techniques.  The  chemical  composition 
of  the  particular  alloys  investigated  (Pbggln.^  an<^  Pbg^Ingg)  WaS 
determined  using  residual  resistance  ratio  methods. 

The  pressure  dependence  of  the  energy  gap  and  transition  temp¬ 
erature  were  studied  and  found  to  be  in  excellent  agreement  with 
the  recent  data  of  Hansen  et  al  (1973)  and  Galkin  et  al  (1973). 

The  superconducting  energy  gap  was  found  to  decrease  at  a  greater 
rate  than  the  transition  temperature  indicating  a  trend  towards 
weaker  coupling.  The  magnitude  of  this  trend  was  determined  and 
found  to  be  in  good  agreement  with  the  theoretical  prediction  of 
Trofimenkoff  and  Carbotte  (1970). 

The  pressure  dependence  of  the  principle  phonon  peaks  of  the 
system  was  also  investigated  including  the  localised  mode  first 
observed  by  Rowell  et  al  (1965)  ,  and  the  corresponding  Gruneisen 
mode  gammas  for  these  peaks  calculated.  Further  information 
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regarding  the  pressure  dependence  of  the  phonon  spectrum  was  obtained 
using  the  McMillan  inversion  programme  (1965)  for  the  Pb^In^  alloy. 
Data  extracted  by  this  process  were  used  in  conjunction  with  the 
McMillan  equation  for  (1968)  from  which  fairly  good  agreement  of 
the  predicted  pressure  dependence  of  with  the  directly  obtained 
experimental  value  was  in  evidence. 

By  combining  existing  experimental  data  on  the  density  of  states 
at  the  Fermi  surface  with  experimental  results  obtained  here  for  the 
re-normalisation  function,  the  pressure  dependence  of  the  band 
structure  density  of  states  has  been  derived.  The  values  so  obtained 
are  in  obvious  disagreement  with  that  of  the  free  electron  gas.  Based 
on  existing  experimental  evidence  (Anderson  and  Gold  (1965)  and 
Anderson  et  al  (1966))  and  simple  assumptions  therefore,  a  value  for 
the  pressure  dependence  of  the  band  structure  density  of  states  at 
the  Fermi  surface  has  been  derived  for  the  case  of  pure  lead  (which 
is  not  expected  to  differ  greatly  from  the  case  at  hand) . 

The  use  of  tunnel  junctions  as  sensitive  pressure  monitoring 
devices  is  strongly  advocated.  Their  response  to  pressure  has  been 
found  to  be  sensitive  and  reproducible  while  being  almost  independent 
of  temperature. 
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CHAPTER  I 


INTRODUCTION 


The  intent  of  the  author  in  writing  this  chapter  is  to  bridge 
the  gap  between  the  Table  of  Contents  and  the  ensuing  thesis 
material . 

An  attempt  has  been  made  in  Chapter  II  to  develop  the  theory 
which  underlies  only  those  aspects  relevant  to  the  understanding  of 
the  experimental  measurements.  The  chapter  begins  therefore  (after 
a  very  brief  historical  introduction  to  superconductivity)  with  an 
account  of  the  theory  leading  up  to  the  development  of  the  super¬ 
conducting  transition  temperature  and  energy  gap  for  which  a  many  body 
approach  utilising  the  Bogoliubov— Valatin  transformations  has  been 
chosen  in  favour  of  the  original  (and  more  usually  employed)  many  body 
theory  of  Bardeen,  Cooper  and  Schrieffer  (1957),  henceforth  B.C.S. 
Frequent  recourse  is  however  made  to  the  B.C.S.  theory.  The  advantage 
of  this  choice  is  that  it  can  be  extended  to  encompass  the  theory  of 
electron  tunnelling  which  provided  the  analytical  tool  for  the  study  at 
hand.  Extension  of  the  theory  to  cover  the  case  of  the  so  called  strong 
coupling  superconductors  follows.  The  principle  difference  in  the 
phonon  spectrum  of  the  alloys  from  that  of  pure  lead  is  the  existence 
of  a  localised  mode  of  vibration.  The  final  two  sections  of  Chapter 
II  consist  firstly  of  an  attempt  to  account  for  the  observation  of 
phonon  processes  through  electron  tunnelling  and  secondly  a  simplistic 
justification  for  the  existence  of  the  localised  mode. 
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Chapters  III  and  IV  contain  descriptions  of  the  experimental 
apparatus  employed,  the  latter  being  devoted  entirely  to  the  electronic 
detection  systems  (and  their  operation).  Circuit  analyses  are 
presented  verifying  the  linear  response  of  the  two  bridge  circuits  that 
were  used. 

Experimental  method  is  the  subject  of  the  fifth  chapter. 

In  Chapter  VI,  with  the  exception  of  section  five,  ordering 
of  the  results  (and  discussion  thereof)  has  been  based  on  the 
experimental  ease  with  which  they  were  obtained  and  more  importantly, 
on  the  dependence  of  subsequent  sections  upon  their  predecessors. 


CHAPTER  II 


THEORETICAL  CONSIDERATIONS 


1 .  Towards  the  Microscopic  Theories  of  B.C.S.  and  Bogoliubov 

Before  developing  in  detail  the  theory  which  underlies  the  measure¬ 
ments  made,  it  is  felt  pertinent  to  very  briefly  acquaint  the  reader 
with  some  of  the  historical  background  leading  up  to  the  successful 
theories  of  Bardeen,  Cooper  and  Schrieffer  (1957)  henceforth  B.C.S. 
and  of  Bogoliubov  (1958)  the  latter  being,  in  the  main,  the  basis  of 
approach  used  thereafter. 

The  phenomenon  of  superconductivity  was  discovered  in  1911  by 
K.  Onnes  whilst  studying  the  electrical  resistance  of  mercury  at  low 
temperatures  when  he  noticed  that  at  some  particular  temperature,  the 
resistance  of  his  specimen  dropped  to  a  value  which  he  could  not 
differentiate  from  zero.  This  temperature  came  to  be  known  as  the 
transition  (or  critical)  temperature  T^  and  is  probably  the  best  known 
property  characteristic  of  the  superconducting  state.  Subsequent 
studies  showed  this  state  to  exist  in  numerous  metals  and  alloys 
and  that  it  was  reversible  in  the  sense  of  other  thermodynamic  phase 
changes.  Attempts  to  explain  the  physical  basis  for  the  existence 
of  this  new  state  were  first  proposed  by  Gorter  and  Casimir  (1934)  and 
by  F.  and  H.  London  (1934-1935).  Both  of  these  models  were  of  a 
phenomenological  nature  treating  the  electrons  in  the  metal  as  two 
fluids,  one  component  being  ascribed  the  properties  of  normal  electrons 
whilst  the  other  gave  rise  to  the  superconducting  properties.  The 
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difference  in  the  approaches  was  that  Gorter  and  Casimir  used  a 
thermodynamical  treatment  whilst  the  London  treatment  was  electro¬ 
dynamic  in  nature. 

The  first  big  step  towards  an  understanding  of  the  superconducting 

state  at  the  microscopic  level  was  achieved  by  H.  Frohlich  (1950)  who 

proposed  that  the  critical  temperature  of  a  superconductor  was  related 

_i 

to  its  isotopic  mass  according  to  T  ^  M  2  and  thus  indicating  that  the 
lattice  was  instrumental  in  the  onset  of  superconductivity.  In  arriving 
at  this  result,  Frohlich  had  considered  the  interaction  of  electrons 
with  the  lattice  (virtual  phonon  exchange)  and  concluded  that  if  such 
interactions,  which  were  attractive,  could  overcome  the  normal 
electrostatic  repulsion  between  these  electrons  then  the  superconducting 
state  would  be  realised.  This  attempt  proved  quantitatively  inadequate 
but  the  lattice  involvement  had  been  arrived  at  which  was  the  key  to  the 
eventual  successful  microscopic  theories. 

In  1956,  Cooper  showed  that  two  electrons  of  opposite  momenta  and 
spin  interacting  via  a  momentum  dependant  potential  formed  a  bound  state 
above  the  Fermi  surface  if  the  interaction  between  them  was  attractive 
no  matter  how  weak.  At  this  time,  evidence  was  growing  suggesting 
the  existence  of  an  energy  gap  in  the  single  particle  electron  density 
of  states  near  the  Fermi  surface. 

Based  mainly  on  the  mechanism  proposed  by  Fr’ohlich  and  the  refine¬ 
ments  of  Cooper,  invoking  the  pair  concept,  an  extensive  and  successful 
theory  of  the  superconducting  state  was  developed  by  Bardeen,  Cooper  and 
Schrieffer  in  1957  and,  using  a  different  mathematical  approach,  by 
Bogoliubov  in  1958  from  which  such  properties  as  the  transition 
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temperature  and  the  suggested  single  particle  energy  gap  naturally 
arose.  These  theories  were  later  extended  to  describe  the  so  called 
strong-coupling  superconductors  which  were  found  to  exhibit  structure 
in  the  electronic  density  of  states  curve  which  could  not  be  accounted 
for  by  either  of  the  above  mentioned  theories.  These  and  other 
deviations  from  the  B.C.S.  predictions  have  been  attributed  to  the 
specific  nature  of  the  lattice  involvement  with  superconductivity  in 
which  for  some  classes  of  superconductors  the  electron  phonon  interaction 
is  weak  and  can  be  described  by  the  B.C.S.  or  Bogoliubov  theories  whilst  for 
others  the  interaction  is  strong  and  necessitates  a  different  approach. 

That  class  of  superconductors  which  can  be  adequately  described  by  the 
theories  of  B.C.S.  or  Bogoliubov  have  therefore  become  known  as  weak 
coupling  (or  B.C.S.)  superconductors  as  opposed  to  those  in  which  the 
lattice  involvement  is  more  prominent  which  are  known  as  strong- 
coupling  superconductors. 

2 .  The  B.C.S.  Hamiltonian 

In  an  attempt  to  arrive  at  a  pervasive  theory  for  superconductivity 
the  search  leant  towards  finding  an  interaction  between  electrons  which 
would  be  attractive  since  this  would  lead  to  a  lowering  of  the  total 
energy  of  the  superconducting  state  with  respect  to  that  of  the  normal 
state  as  was  being  found  experimentally.  In  view  of  the  growing  belief 
that  the  lattice  was  involved  with  the  onset  of  superconductivity  and 
with  the  success  of  the  isotope  effect,  the  Hamiltonian  of  Frohlich 
(from  which  the  isotope  effect  evolved)  seemed  an  obvious  starting 
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point  for  the  development  of  the  theory.  This  Hamiltonian  may  be 


written  explicitly  in  the  form 


2-1 


where  W, 


,  ,  is  a  matrix  element  of  the  form 
kk  q 


(ek  -  W  -  (%)2 


2-2 


represents  the  system  of  noninteracting  electrons  and  phonons  whilst 


the  succeeding  term  represents  the  interaction  of  electrons  with  the 
lattice.  In  physical  terms  we  may  think  of  the  interactive  term  as 
representing  a  self-energy  associated  with  the  electrons  in  the 
following  way.  An  electron  moving  through  the  lattice  distorts  the 
lattice  (ion  sites),  the  lattice  in  turn  then  reacts  on  the  electron 
due  to  its  changed  electrostatic  field.  As  the  electron  passes  by, 
the  lattice  is  allowed  to  relax.  Such  distortions  of  the  lattice  are 
quantised  in  terms  of  phonons  and  we  may  therefore  think  of  the 
electrons  moving  through  the  lattice  constantly  absorbing  and  re-emitting 
phonons.  This  surrounding  of  the  electrons  by  phonons  has  given  rise 
to  the  concept  of  "dressed"  or  "clothed"  electrons.  The  ability  of  such 
a  dressed  electron  to  repel  a  second  electron  in  its  vicinity  is 
therefore  reduced  by  this  positive  screening  effect  and  it  is  this 
reduction  of  its  potential  to  repel  that  is  referred  to  as  the  self¬ 
energy  of  the  electron.  It  should  be  noted  that  the  exchanged  phonons 
are  very  short  lived  and  that  as  a  consequence  of  the  uncertainty 
principle,  energy  need  not  be  conserved  in  the  process,  to  this  end 
the  phonons  involved  are  virtual. 


Frohlich's  proposal  therefore  was  that  the  desired  attractive 
interaction  could  be  realised  if  this  screening  of  the  conduction 
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electrons  by  the  positive  ion  lattice  sites  was  such  as  to  invoke  an 
effective  attractive  force  between  electrons  which  was  larger  than  their 
ability  to  repel  each  other  as  bare  electrons.  The  lowering  in  total 
energy  (condensation  energy)  associated  with  the  superconducting 
phase  change  as  predicted  by  this  model  turned  out  to  be  three  or  four 
orders  of  magnitude  greater  than  that  observed.  It  can  be  seen  from  2-1 
that  the  Frohlich  Hamiltonian  included  interactions  between  electrons  of 
all  momenta  and  either  spin  and  was  in  fact  the  reason  for  the  over 
estimate  of  the  condensation  energy. 

A  refinement  proposed  by  Cooper  in  1956  eventually  led  to  an 
acceptable  form  for  the  Hamiltonian  when  he  showed  that  if  there  was  a 
net  attraction  between  a  pair  or  electrons  just  above  the  Fermi  surface 
then  they  could  form  a  bound  state.  He  also  showed  that  such  electrons 
which  could  form  this  bound  state  by  phonon  mediation  must  be  within 
an  energy  range  of  the  Fermi  surface  where  is  of  the  order 

of  an  average  phonon  energy  of  the  metal,  and  that  for  maximum  lowering 
of  energy  all  such  pairs  must  have  the  same  total  momentum.  This 
model  may  be  visualised  as  follows.  One  electron  moving  through  the 
lattice  distorts  the  lattice  thereby  emitting  a  phonon.  The  second 
electron  in  the  vicinity  is  then  affected  by  this  disturbance  (phonon 
absorption)  as  it  propagates  through  the  lattice.  If  the  energy  of 
the  phonon  which  mediates  the  process  is  greater  than  the  electronic 
energy  change  invoked  thereby,  then  the  resulting  interaction  is 
attractive.  Such  pairs  of  electrons  are  referred  to  as  Cooper  pairs. 

This  pair  concept  was  extended  to  the  macroscopic  occupation 
level  by  Bardeen,  Cooper  and  Schrieffer  who  formulated  a  reduced 
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problem  in  which  only  pairs  whose  constituents  had  equal  and  opposite 
momentum  were  considered.  This  reduces  the  Frohlich  Hamiltonian  to:- 

2-3 


HBCS  ^£kCksCks  2  Vkk' Ck' s '  C-k' sC-ksCks’ 
ks  kk  s 


where  \k'  =  -2w-k,k,k'-k‘ukk' 


2-4 


U,  ,  is  a  screened  Coulomb  repulsion  term  added  to  the  Frohlich 
K.K. 

Hamiltonian  of  2-1.  If  V^t  is  positive  then  the  electron-electron 
interaction  is  attractive. 

It  can  be  shown  that  the  ground  state  energy  described  by  this 

Hamiltonian  is  minimised  if,  in  addition  to  pairing  only  electrons  whose 

momenta  are  equal  and  opposite,  we  also  require  that  such  paired  electrons 

have  opposite  spin.  Thus  by  imposing  this  restriction  and  applying 

the  convention  that  an  operator  preceded  by  a  minus  sign  implies  a 

spin  down  state  whilst  one  with  a  positive  sign  implies  a  spin  up 
4*  4* 

state  e.g.  c^,  =  c^tf  ;  c_^f  =  c(-k')l’  We  Can  wr^te  t^rie  Hamiltonian 


of  B.C.S.  as 


hbcs  =  ^ek(ckck  +  C-kC-k)  '  ,  J,vkk'ck'c-k'c-kck 
k  kk' 


2-5 


It  is  at  this  juncture  that  we  depart  from  the  B.C.S.  development 
and  proceed  instead  to  study  the  eigenvalues  and  eigenstates  of  this 
Hamiltonian  according  to  the  Bogoliubov  method. 


3.  The  Bogoliubov  Valatin  Formalism 

We  next  proceed  to  diagonalise  the  Hamiltonian  given  by  2-5  by 
first  introducing  a  transformation  to  a  new  operator  system  defined 


by 


■ 
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t 

Yk  = 

and  their  conjugates 

f  t 

Yk  =  ukck  -  vkc_k 


Y_k  -  V-k  +  Vk 


t 


.t 


Y-k  =  ukC-k  +  Vk 


2-6 


u  and  v  are  chosen  to  be  real  constants  related  by  the  condition 


V  +  V  - 1 


2-7 


Such  a  choice  ensures  that  the  newly  defined  operators  satisfy  the 


Fermi  anticommutation  laws  namely :- 

[AT,A+]^  =  [A  ,A  ]  =  0 

L  p*  qJ+  L  p  qJ  + 


[A  ,A+K  =  6 

1  p  qJ+  pq 


2-8 


where  the  A  represents  either  the  newly  created  operators  or  the 
traditional  electron  (c  type)  annihilation  and  creation  operators. 

In  order  to  justify  the  relationship  2-7  one  such  commutation  relation 
for  the  new  operator  system  is  shown  here. 

IV*k3  =  (ukck'vkc-k)  (ukcI’vkc-k)  +  <  vl'V-k) ( Vk-ykctk) 

multiplying  out  and  regrouping  terms  yields 

=  uk2[ck,ck]+  -  ukvk[Vc_k]  +  ~  ukvk[  Jk , ck]  +  +  vk2  [c^k,c_kJ  + 


V  +  V  E  1 


using  assumption  2-7 
The  inverses  of  the  B.V.  transformations  are  easily  obtained  in  the 


following  form:- 

\ 

ck  =  Vk  +  V-k  cl  =  V+k  +  v-k 

2-9 

c-k  '  v-k  ■  Vk  c-k  ‘  v-k  ■  Vk 

We  next  proceed  to  transform  the  B.C.S.  Hamiltonian  to  our  new  operator 
system.  Writing  2-5  as  the  sum  of  kinetic  and  interaction  terms  we 


have :  - 


■ 
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H  =  Ht  +  Hy  =  (in  the  transformed  system) 


2-10 


Considering  first  the  term  H^,  we  have:- 

+  + 

Ht  =  +  c_^c_^) 

k 

=  E£k{  ( vit  +  vkY-k>  (uk\  +  V-k>  +  ( V-k  -  Vk> 

k 

x(ukY-k  '  vl> } 

=  ,E£k[uk2\\  +  VW-k  +  Vk(Yk^-k  +Y-kV  +  vkVk  + 

k 

Vw-k  -  +  \Y-k)J 


2-11 


For  simplicity  we  write  this  expression  term  by  term  as 

H  E  Ze  [P  +  Q  +  R+  S  +  T  +  U]  then  using  the  anticommutation  laws 

T  k  k 


for  the  y 


t  t 


R  +  U  =  2ukvk(ykY_k  +  Y_kYk) 

Q  -  VY.ak  -  v,  2(1  -  yIy.J 


\  J-k'-k  vk  VJ_  '-k'-k2 

s  =  Vvl  =  vk2(1  -  YkV 


p  +  s  ■  V  +  (uk2  -  vk2)\Yk 

I  +  Q  =  vk2  +  (uk2  -  vk2)Y!kY.k 


before  re-collecting  these  terms,  new  number  operators  are  defined 
for  the  gammas  in  analogy  with  the  electron  number  operators  i.e. 

Yk\  E  "V  and  Y-kY-k  =  m-k 

We  can  thus  write  HT  in  the  following  way, 


ht  =  Sek[(p  +  S  +  T  +  Q)  +  (R  +  U)J 

K 


to  get:- 


11 


HT  =  ?£k[2vk2  +  (uk2~vk2)(mk  +  m-k)  +  2ukvk(\>tk+Y_kYk>]  2-12 
k 

The  interaction  term  H^  of  equation  2-10  is  treated  in  a  similar  fashion 
(more  extensive  use  is  made  of  the  commutation  relations  and  the 
grouping  of  terms  requires  greater  insight  due  to  the  complexity  of 
this  term)  to  yield 

^  =  -  hvkk’  [uk'vk’ukvk(1  -  v  ~  m-k,)(1  ■  "v  ■  m-k)+ 

kk 

VV(1  -  V  -  m-k’)(uk2  -  \!)(Vk  +  V-k)J  + 

(fourth  order  off  diagonal  terms)  2-13 

The  Hamiltonian  is  now  diagonalised  for  the  ground  state  of  the 
superconductor  for  which  we  assume  that  all  the  occupation  numbers 
such  as  m.  vanish.  This  is  achieved  by  requiring  that  the  sum  of 
off  diagonal  contributions  to  H  ^  from  and  cancel  one  another 


l.  e. 


pwk(v-k +  Y-kV 

k 


3,vkk-w(uk2-  vk2)(v-k +  Vi’ 

kk 


=  0 


2-14 


in  which  the  assumption  is  made  that  the  expectation  value  of  the 
fourth  order  off  diagonal  elements  is  negligible  in  the  ground  state 
of  the  system.  Elimination  of  the  off  diagonal  elements  in  this  manner 
enables  us  to  consider  the  system  as  being  composed  of  independent 
fermions.  Considering  equation  2-14  term  by  term  we  get:- 


2ekukvk  -  (uk2  -  vk2)^,\k'uk'vk'  =  0 

k 


2-15 


which  dictates  the  form  that  u  and  v  (the  B.V.  transformation)  must 
have  to  be  consistent  with  the  diagonalisation  condition.  Since 
u  and  v  are  interrelated  they  may  be  defined  in  terms  of  a  single 
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parameter  as  follows :- 


and 


v. 


=  (h  + 


V' 


2-16 


Equation  2-15  then  becomes 
2ek^  "  Xk2')  2  +  2Xk  Ak 

where  we  have  introduced  a  new  quantity  defined  by:- 

\  ■  IM1  -  v2)% 

k 


Thus  x,  = 
k 


±e. 


2(£kz+  V)J* 


substitution  for  in  2-18  yields  a  form  for  A^:- 


\  ■  ^,\k- 
k 


V 


ten,.*  +  A  J)% 


2-17 


2-18 

2-19 


2-20 


This  latter  equation  represents  an  integral  equation  which  can  in 

principle  be  solved  if  we  know  the  form  of  the  interaction  potential 

V  . .  In  order  to  conserve  the  total  number  of  particles  in  the 
kk 

system  we  must  choose  the  Fermi  surface  as  the  zero  level  of  energy 
for  the  system.  The  proof  of  this  is  seen  below  which  also  enables 
us  to  decide  upon  the  sign  in  equation  2-19.  The  number  of  electrons 
in  the  system  is  given  by:- 

N  =  E(ckCk  +  c-kC-k^ 
k 


The  expectation  value  of  N  is,  using  equation  2-12:— 

<o | N | o>  =  N<o | o>  =  <N>  =  £<o | 2v^2 | o>  +  <o | (uk2  -  vk2) 

k 

x(mk  +  m-k)l°>  +  <0l2ukvk(\Y-k  +  Y-kYk)l°>  2' 

Now  for  the  ground  state  of  the  system  the  =  0  in  the  second  term 


- 
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whilst  the  third  term  can  be  shown  to  produce  only  zero  eigenvalues 
by  replacing  the  y  according  to  equation  2-6.  Thus  the  expectation 
value  of  N  reduces  to:- 

<N>  =  E2vk2  =  2(1  +  2xk)  2-22 

k  k 

If  we  had  no  interactions  i.e.  all  the  electrons  resided  below  the 
Fermi  surface  then  we  know  that  N  must  have  the  form 


<N>  =  £  2 

k<k„ 


and 


<N>  =  0  for  k>k. 


2-23 


where  the  factor  of  two  arises  because  of  the  dual  spin  nature  of 
fermions.  For  the  non  interacting  case  then  we  must  put  V^,  =  0 
and  thus  from  equation  2-18,  =  0  also.  By  substitution  then  in 

equation  2-19,  we  see  for  the  non  interacting  case  that:- 
xk  =  ±% 

Where  the  positive  sign  is  chosen  to  represent  the  system  below  the 
Fermi  surface  and  the  negative  sign  for  above  the  Fermi  surface 
thereby  making  equation  2-22  consistent  with  the  real  physical 
situation  of  equation  2-23.  So  for  the  normal  metal  the  B.V. 


transformation  has  the  form 

v, 


x. 


E, 


Where  e,  here  and  henceforth  will  be  used  to  denote  measurements 
k 

with  respect  to  the  Fermi  surface.  Turning  attention  now  to  the 


' 
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interacting  case  (V^,^  0)  of  a  superconductor  we  can  see  from 
equation  2-22  that  for  the  number  of  particles  to  be  conserved 
that  must  be  an  odd  function  of  (e^  -  e^)  as  was  the  case  above 
and  thus  we  choose  the  negative  sign  for  and  have  therefore :- 

x,  =  2-24 

+  V* 


Notice  that  V^,  ->0  (and  hence  -*  0)  has  the  highly  desirous  effect 

of  encompassing  the  normal  metal  case.  Equation  2-24  describes  the 

form  that  the  B.V.  transf ormatiors (i.e.  the  restriction  on  our 

original  u  and  v)  must  have  for  the  Hamiltonian  to  be  of  diagonal 

form  and  are  shown  pictorially  below. 

v  x. 


f  k 


Bardeen,  Cooper  and  Schrieffer  solved  equation  2-20  for  delta  using 

a  very  simple  potential  for  which  they  proposed  the  form 

Vkk' 

a 


V. 


kk* 


V  for  |e .  |  &  |ek. 


0  otherwise 


~1 

< — TiW  — » 

L 

D 

->  e, 


Expressed  in  integral  form,  equation  2-20  becomes 


-  u 


D(ek,) 


V 


(ek«2+  >?s 


Vkk’dek’ 


2-25 


— OO 
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Since  V,  ,  .  is  constant  and  so  is  Tioo_,  (the  Debeye  energy)  then  it 
kk  D 

follows  that  A,  is  constant  and  we  may  therefore  omit  the  subscripts 

K. 

and  take  V  outside  the  integral.  D(e^,)  represents  the  density  of 
electron  states  (for  electrons  of  one  spin)  and  is  assumed  constant 
at  the  Fermi  surface  thus  equation  2-25  reduces  to 

r% 


=  VD 
2 


A 


de 


(e2  +  Az) 


2\h 


-Iko. 


D 


This  integral  can  now  be  easily  solved  using  the  substitution 
e  =  A  sinh  x  to  yield 

A  =  'ha)D/{sinh(l/DV)  }  2-2' 

For  a  great  many  metals,  estimates  of  the  quantity  (DV)  yield  values 
^  10  ^  whence  it  is  easy  to  see  from  the  definition  of  the  hyperbolic 
function 

sinh  x  =  %(e  -  e  ) 

that  we  may  reasonably  make  an  approximation  to  2-26  of 

A  =  2hwDexp[-l/VD(ef)]  2-2 

Those  metals  for  which  the  above  approximation  is  valid  are  said  to 
be  weak-coupling  superconductors.  It  can  be  seen  from  2-27  that  A 
is  an  extremely  small  quantity  and  that  because  of  this  the  parameter 
x  defined  by  equation  2-24  differs  from  ±h  only  very  close 

K. 

(within  ±A)  of  the  Fermi  surface.  Examining  this  implication  from 
a  wider  viewpoint  we  see  that  as  a  consequence  of  x^  having  the 
values  in  regions  remote  from  the  Fermi  surface,  the  values  of 
u  and  v  become  either  unity  or  zero  and  hence  from  equation  2-6 
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we  see  that  our  newly  introduced  operators  gamma  revert  back  to  a 
single  c  type  operator  as  in  the  case  of  a  normal  metal. 

We  have  to  this  point  then  built  up  a  picture  of  a  superconductor 
which  resembles  a  normal  metal  except  within  a  small  region  bounding 
the  Fermi  surface  where  the  electrons  interact  in  a  very  specific 
manner.  The  nature  of  this  interaction  is  the  pair  concept  alluded 
to  earlier  and  the  gamma  operators  set  forth  at  the  beginning  of 
this  section  are  those  which  create  and  destroy  such  Cooper  pairs. 

The  term  quasiparticle  is  used  to  describe  a  single  electron  while 
in  a  paired  state.  The  parameter  A  is  associated  with  the 
energy  of  formation  (or  destruction)  of  such  pairs.  It  is  the  very 
small  lowering  of  energy  associated  with  the  pairing  of  electrons 
combined  with  the  fact  that  this  pairing  only  occurs  in  a  very 
restricted  region  at  the  Fermi  surface  that  accounts  so  nicely  for 
the  experimentally  observed  small  lowering  of  the  total  energy  of 
the  metal  when  in  the  superconducting  state. 


4 .  The  Superconducting  Ground  State  and  Energy  Gap  (T  0) 


Having  established  in  the  previous  section  the  condition  that 
the  off  diagonal  elements  of  HBV  =  +  Hy  vanish,  we  are  left  with  the 

diagonalised  Hamiltonian 

hbv  =  Eek[2vk2  +(uk2-  vk2)(mk  +  m-k)J"E.  ,vkk'Wukvk 

k 


u  -  v  -  m-k')(1  -  “k  -  m-k) 


2-28 


which  in  the  ground  state  where  the  m  =  0  reduces  to 
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2-29 


The  ground  state  wave  function  ¥  for  the  superconductor  must  be  an 
eigenfunction  of  the  diagonalised  Hamiltonian  so  that  for  all 
possible  k:- 

v,  k  >  =  Y1|'i'>  =  0  2-. 

'k1  o  -k1  o 

This  in  effect  is  an  attempt  to  create  a  state  lower  than  the  ground 
state  defined  by  ¥  .  In  order  to  find  the  form  of  T  the  following 
argument  is  extended.  Since  the  gammas  obey  the  usual  Fermi  anti¬ 
commutation  laws,  products  such  as  =  Y.^Y.^  =  0  ancl  so  equation 

2-30  is  satisfied  by  operating  on  the  vacuum  state  with  all  of  the 
possible  products  of  the  operators  yk  and  because  such  products 

would  always  contain  a  product  of  the  type  ykYk  and  thereby  yield 
a  zero  eigenvalue.  For  all  possible  k  values  therefore  this  may  be 
stated  as 


2-31 


where  again  use  has  been  made  of  the  original  definition  of  the 
gammas.  This  expression  simplifies  by  using  the  anticommutation 
relations  and  normalising  throughout  by  the  product  of  all  of  the 


v,  to  give  for  the  ground  state  wave  function 

K. 


2-32 


k 


It  is  possible  to  justify  the  above  expression  on  a  purely  mathematical 
basis  in  which  we  first  propose  the  form  for  \VQ>  as  given  in  2-32 
and  show  that  it  is  consistent  with  the  physical  proposal  of  equation 
2-30,  the  derivation  for  which  will  not  be  presented  here. 
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Quasiparticle  excitations  of  the  system  can  be  created  simply 
by  acting  upon  the  ground  state  wave  function  with  the  creation 
operators  and 

The  ground  state  energy  of  the  superconducting  state  £g  can  be 
obtained  by  calculating  the  expectation  value  of  the  Hamiltonian 
given  by  equation  2-29. 

e  =  <K„>  =  <o  I  E2£.  v,  2  -  E  V,  i  , u.  ,  v,  , u.  v  j  o> 


BVG 


'kk 


kk 


, vkk,lV  vk’~kvk' 


Using  the  model  potential  proposed  by  B.C.S.  we  see  that  all  of  the 
quantities  in  the  above  expression  are  constants  and  therefore 


=  E2ekV  "  E-V-uv'vk'u 


V. 


kk 


, vkk,uk' vk,ukvk 


2-33 


When  excitations  are  present  however  (m  £  0)  we  must  revert  to  the 
Hamiltonian  of  equation  2-28  the  terms  of  which  may  be  rearranged 
to  give 


£ 

s 


+ 


E  (iy  +  m_k)  [  (uk2 


V)Sk 


+  2 


Yk'uk-vk'J 


+  higher  order  terms 


2-34 


and  where  use  has  been  made  of  the  expression  derived  for  £g. 

By  substituting  in  equation  2-34  (using  equations  2-16 ,  2-20  and 
this  expression  can  be  rewritten  in  the  form:- 

HBV  =  es  +  E(V+  A2>"2(mk  +  ”-k)  + . 


2-22) 


2-35 


It  can  now  be  seen  that  this  represents  a  system  whose  ground  state 

energy  is  given  by  £  and  whose  (quasiparticle)  excitations  are  given 

s 

by:- 

\  -  <V  +  a2>% 

These  quasiparticles  cannot  be  created  singly  since  any  physical 


2-36 
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process  is  represented  by  at  least  two  electron  operators  e.g. 
scattering,  which  results  in  two  gamma  type  operators.  The  minimum 
energy  required  to  create  (or  break  up)  such  a  pair  of  quasiparticles 
is  2A  as  can  be  seen  from  equation  2-36.  When  a  Cooper  pair  is 
broken  up  one  electron  can  be  excited  across  the  Fermi  surface  to  a 
state  (e  +  A)  whilst  the  other  falls  back  into  the  Fermi  sea  (ef  -  A) 
until  re-pairing  with  other  electrons  occurs.  In  this  way  the  energy 
gap  in  the  single  particle  electron  density  of  states  that  is  observed 

can  be  visualised. 

In  analogy  with  the  density  of  states  of  a  normal  metal,  the 
density  of  quasiparticle  states  for  the  superconductor  can  be 


deduced  from 


D (E)  =  dk 
dE 


d|  dk 
dE  d£ 


where  dk/d£  is  the  density  of  Bloch  states  of  a  normal  metal, 
Using  equation  2-36  we  find  that:- 


D(E)  = 


2D(ef)_ 


for  |E|>A 


(E2  -  A2) 


2\h 


0  for  |e  I^A 

Pictorially  this  compares  with  the  normal  metal  as  below 


Normal 


2-37 
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E 
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A 


/ 
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/ 


/ 


/ 


Superconducting 


5 .  Finite  Temperature  Considerations  and  the  Transition  Temperature 

In  the  previous  section  it  had  been  assumed  that  the  number  of 
quasiparticle  excitations was  very  small,  i.e.  A»kT,  a  situation  we 
might  envisage  as  T  ->  0.  We  now  treat  the  more  physical  case  of 
T  ^  0  where  we  no  longer  assume  the  number  of  quasiparticle 
excitations  to  be  negligible.  In  so  doing  we  will  arrive  at  a 
modified  form  for  the  quasiparticle  excitation  energies  and  energy  gap,  and 
derive  an  expression  for  the  transition  temperature.  In 
proceeding,  we  assume  that  the  Hamiltonian  for  the  system  can  still 
be  diagonalised  (an  assumption  to  be  justified).  We  are  left  then 
with  the  Hamiltonian  given  by  equation  2-28.  The  energy  required  to 
create  a  quasiparticle  excitation  will  be  given  by:- 


2-38 


If  we  apply  this  to  equation  2-28  we  get 


!  V,  A  vii  I  Ui  |  V, 

k  L  f  kk  k  k 


In  defining  the  excitation  energies  as  in  2-38  we  may  treat  the  problem 
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as  a  system  of  independent  fermions  and  obtain  an  approximate 

solution  by  replacing  the  m  by  their  thermal  averages  according 

to  Fermi-Dirac  statistics. 

<hl  >  =  <m_,  >  =  _ 1 _  =  f  (Ei) 

exp(E^/kT)  4-  1 


The  off  diagonal  terms  are  treated  in  exactly  the  same  way  as  in  the 
previous  case  (T  -*■  0)  except  that  we  do  not  now  assume  that  the  m 
vanish  also,  but  rather  they  take  the  form  of  their  thermal  average 
as  in  equation  2-40.  The  resulting  expression  is:- 

2-41 


2£kVk  ■  (uk2  -  vk2)^,vkk'uk'vk,[1  “  2f(V]  0 


Comparing  this  with  equation  2—15  we  see  that  they  differ  only  by  the 
factor  in  the  square  brackets  containing  the  Fermi  function. 


We  define  a  modified  gap  function  here  which  is  temperature  dependent. 


VT)  ■  J,vkk'ft-  v 

k 


2[1  -  2f(Ek,)] 
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The  development  is  now  as  for  the  case  of  T  +  0  (with  the  B.V. 
transformation  x  being  temperature  dependent)  and  leads  to  the  integral 
gap  equation  for  non  zero  T  of:- 


A,  (T)  =  JgZ  V 


V 


v<T> 


[1  -  2f(Efc,)] 


[A 

e 


k' 


+  A 


k' 


(T)] 
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By  substituting  into  the  expression  for  (equation  2-39)  as  before 
and  including  the  explicit  form  of  the  Fermi  functions  given  by 
equation  (2—40)  we  get  for  the  quasiparticle  energies 

Ek  =  [ek2  +  A2(T)]^  2-44 

This  expression  is  exactly  that  obtained  earlier  but  with  the  temperature 


dependent  gap  function. 
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Using  the  model  potential  of  B.C.S.  and  solving  as  before 


the  integral  equation  2-43  for  the  gap  we  find  that 

f  1 

S~2 


VD(ef) 


0 


D 


tanh[(P  +  A2)^/2kT]d£  =  1 


2-43 


(£2  +  A2) 


2  \"2 


This  clearly  is  a  very-  complicated  implicit  expression  for  the 
energy  gap  at  finite  temperature.  It  can  be  seen  however  that  as 
T  increases,  the  numerator  decreases  and  therefore  for  the  left 
hand  side  of  this  expression  to  remain  constant  equal  to  unity, 
the  denominator  must  decrease  also.  This  implies  therefore  that  A 
is  a  monotonically  decreasing  function  of  temperature.  Such  is 
the  case  observed  experimentally  as  shown  below;  typically 


It  is  very  satisfying  that  if  T  -  0  then  this  expression  yields 
identically  the  previous  expression  derived. 

To  obtain  an  expression  for  the  transition  temperature  we 
need  to  impose  the  condition  that  A  -*■  0  as  T  -►  Tc  in  this  integral 
which  can  then  be  performed  to  yield 

kT^  =  1*  14fi(jo^exp  [-l/VD(e^.)  J  2-46 

having  made  use  of  the  weak  coupling  approximation  kTc«ftwD. 


■ 
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Recalling  the  expression  for  the  T  =  0  case  for  A(0)  (equation  2-27) 

and  comparing  with  this  we  have  that,  for  the  simple  constant  potential 

model  of  B.C.S.  in  the  weak  coupling  assumption :- 

2A(0)  =  3.53  UNIVERSALLY  2-47 

kT 

c 

This  key  result  is  often  used  experimentally  as  an  expression  of  the 
weak  coupling  approximation. 

6 .  Theory  of  Tunnelling  and  the  Density  of  States 

Consideration  is  now  given  to  the  transport  of  electrons  across 
an  insulating  barrier.  This  might  be  represented  by  the  diagram 
below  for  the  case  when  both  metals  are  in  the  normal  state, 


e  e 


where  the  Fermi  level  energy  of  metal  one  has  been  raised  by  an 
amount  eV  with  respect  to  that  of  metal  two.  The  quantities  D^(e) 
and  D  (e)  represent  the  density  of  electron  states  in  metals  one 
and  two  respectively.  Notice  that  in  accordance  with  the  diagram 
on  page  19  that  the  dotted  line  represents  the  density  of  available 


states  in  a  normal  metal  whilst  the  shaded  region  represents  the 
density  of  occupied  states  as  obtained  by  multiplying  the  D(e)  by 
the  respective  Fermi  function. 
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A  Hamiltonian  representing  this  system  may  be  written  as  the 
sum  of  three  terms  the  first  of  which  describes  those  electrons 
of  metal  one  on  the  left  hand  side  and  the  second  for  those  on 
the  right  whilst  the  third  is  a  term  relating  to  the  transfer  of 
electrons  from  side  to  side.  Thus:- 

H  ■  ^kckscks  +  tEekdksdks  +  *  Tkk'  <c+ksdk's  +  dk'scks>  2-48 

ks  ks  kk  s 

Since  it  is  proposed  here  to  calculate  the  net  tunnel  current  we 
need  only  concern  ourselves  with  the  third  term  of  this  expression. 

The  transition  probability  in  lowest  order  is  given  by 

It  2  I  =  I  T I  2 
I  kk'  1  1  1  ’ 

Because  of  the  applied  potential  difference  V,  the  energy 
of  electrons  of  metal  one  are  raised  in  energy  by  an  amount  eV. 

These  electrons  of  energy  can  then  be  elastically  scattered 
(tunnel)  to  an  available  state  k'  in  metal  two  whose  energy  is 
given  by  (e  +  eV) ,  in  this  way  energy  is  conserved.  A  similar 

K. 

argument  applies  to  tunnelling  in  the  opposite  direction.  The 
net  flow  of  electrons  from  left  to  right  is  therefore  proportional 


to 

E  I^UV1  -  -  fk'  (1  -  fk)]S(iv  -  ®k  - eV) 

kk*  s 

where  the  f,  represent  the  probability  of  occupancy  of  a  given 

K. 

state  k.  Expression  2-49  may  be  interpreted  physically  if  we 
associate  f,  with  the  probability  that  the  state  k  from  which  an 

K. 


2-49 


, 
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electron  originates  is  occupied  and  (1  -  fk,)  with  the  probability 
that  the  state  k'  to  which  it  tunnels  is  unoccupied,  energy  conservation 
being  accounted  for  by  the  delta  function.  The  tunnel  current 
therefore  may  be  obtained  from  expression  2-49  by  putting  it  in 
integral  form:- 


INN0C/T2{f  (ek)  "  f(\  +  eV)  }D(ek)D(ek  +  eV)dek 


2-50 


In  obtaining  this  we  have  complied  explicitly  with  the  energy 
conservation  aspect.  T2  is  taken  to  be  the  square  of  some  average 


tunnelling  matrix  element.  If  we  were  to  proceed  from  this  point 


2  2 

using  the  phenomenological  approach  we  would  assume  that  T  and  D 


(arising  from  the  density  of  states  product)  were  constant  and  extricate 
them  from  the  integral,  here  however  a  brief  explanation  is  given  to 
justify  this  procedure.  Using  the  particle  in  a  box  concept  of 
tunnelling  across  a  narrow  potential  barrier  it  can  be  shown  fairly 
easily  that  the  element  T  is  given  by:- 

T  =  Ak  +  B  2-51 

A  and  B  being  constants  and  k  being  the  wave  vector  in  the 
direction  that  tunnelling  takes  place.  The  kinetic  energy 
associated  with  the  state  k  is  simply 


£  =  li2k2  -  £f=>de  =  "ft 2k 


2-52 


Thus  we  see  from  2-51  and  2-52 


T  «  d£ 
dk 


2-53 


But  D (e)  =  dN  =  dN  dk 
d£  dk  d£ 


2-54 


and  we  can  see  therefore  that  the  energy  dependence  of  the  product 
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T  xD  in  the  integral  disappears.  If  we  make  the  justifiable 
assumption  that  dN/dk  is  constant  in  the  region  of  the  Fermi 
surface  we  have  for  I  that 

1  (Ef  )Dr  <£f  )T2/{f  (ek>-f  (6k+eV)  }d6k  2'! 

This  model  is  consistent  with  the  premise  that  tunnelling  takes 
place  at  or  near  the  Fermi  surface  and  that  eV* 

The  factor  f (eL  +  eV)  in  equation  2-55  can  now  be  expanded  as 

K. 

a  Taylor  series  about  whereupon  after  neglecting  terms  containing 
powers  of  (eV)  squared  and  greater  we  have:- 

INN“Dl(£f)Dr(£f)T2/eVf’(ek)dek 

At  low  temperatures,  f  approximates  to  a  delta  function  giving:- 

INN"V£f)Dr(£f)T2eV  2~ 

For  N  -  N  (normal  metal  -  normal  metal)  tunnelling  thus,  an  ohmic 

dependence  is  predicted. 

Attention  is  now  turned  to  the  case  where  one  of  the  metals  is 
in  the  superconducting  state  as  depicted  below. 


£  E 
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Once  again,  in  an  effort  to  derive  an  expression  for  the 
tunnel  current,  we  need  consider  only  the  tunnelling  term  of  the 
Hamiltonian  of  equation  2-48  except  that  those  operators  previously 
assigned  as  d  for  states  on  the  right  hand  side  (metal  two)  must 
now  be  replaced  by  their  quasiparticle  counterparts  in  keeping 
with  the  B.V.  formalism  adopted  earlier.  Explicitly  we  have  then:- 

ht  =  ,?,T[ck(uk'V  +  vk'Y-k')  +  c-k(VY-k'  ■  vk'Yk’) 

kk 


+  <VV  +  VY_ki>ck  +  (uk,y!k.  -  VV)c-k]  2-57 

Notice  that  the  spin  has  been  included  explicitly  here  in  keeping  with 
the  pair  concept  and  has  therefore  been  omitted  from  the  summation 
sign,  it  is  also  very  relevant  to  recall  that  if  the  interaction 
potential  V  ,  responsible  for  the  pairing  tends  to  zero  then  the  y's 
revert  to  single  electron  operators  and  we  have  identically  the 
preceeding  case  for  N  -  N  tunnelling. 

The  physical  interpretation  of  equation  2-57  is  as  follows.  The 
first  term  creates  an  electron  on  the  left  hand  side  of  the  barrier 
and  either  creates  or  destroys  a  quasiparticle  excitation  on  the 
right  (above  or  below  the  Fermi  surface).  The  other  terms  may  be 
interpreted  in  an  analogous  manner. 

JL 

For  a  term  such  as  c^y^,  to  cause  a  real  scattering  (tunnelling) 

we  must  have  that  E^,  =  +  eV  for  an  applied  potential  V.  Note 

that  in  this  case  we  are  measuring  energies  with  respect  to  the 

Fermi  surface  as  was  the  case  in  our  previous  application  of  the 

j*  t 

B.V.  formalism.  Similarly  for  the  term  cky_k, 
scattering  process  we  must  have  -e,  -  eV  =  E,  ,  . 


to  cause  a  real 


' 
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Extending  this  type  of  argument  to  the  remaining  terms  in  2-57 


and  proceeding  as  with  the  N  -  N  case  we  get:- 

5«s‘  £,M2<v2{f(gk)[1  -  f(\')]  ‘  £(Ek')l1  -  f(gk)]}6(Ek'  - 

kk 

+  v,,,2{£(e,)t(E, .,)  -  [1  -  f(e.)][l  -  f(E.  ,)]}«<£.  +  E,  ,  +  eV)) 


2-58 


having  grouped  together  terms  in  u2  and  v2.  By  multiplying  out, 
collecting  terms  and  expressing  in  integral  form  we  have 
INS-/T2D(ek)D(E  =  +  eV)  {uj  [f  (fij-f  +  eV)J  + 


+  vk,2[f(gk)  -  1  +  f(-ek  -  eV)}d£ 


k'  *wk 

k 
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where  the  f  functions  have  been  included  explicitly.  Since  the 
metal  on  the  right  hand  side  of  the  barrier  is  a  superconductor, 
the  density  of  states  function  D(E)  has  the  form  given  by  equation 
2-37.  This  expression  can  be  further  reduced  by  noticing  that 
f (-E)  =  1  -  f (E)  to  give:- 

INScc/T2D(^k^D^  +  eV){f  (ek)  Iuk,  2  +  vk,2J  -  f(ek  +  eV) 


x  [V2  +  vk'2]}dgk 


2-60 


According  to  the  original  B.V.  assumption  u2  +  v2  =  1  hence  we  have 
(replacing  D(e  +  eV)  using  2-37)  that:- 


I  E  |  [f  (€,.)-  f(e,.  +  eV)]d£,. 
<E2  -  A k  k  k 


2-61 


where  the  T2D2  term  again  can  be  taken  outside  the  integral.  This 
is  the  first  point  at  which  the  rigorous  physically  based  theory 
coincides  with  the  phenomenological  approach  most  often  used. 


eV) 


By  noting  that:- 
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f(gk}  “  f(§k  +  eV)  = 


1  for  0  <E  <eV 
0  E>eV  or  E<0 


equation  2-61  becomes  eV 


INs“T2Dl(£f)Dr(£f) 


0 


1  E I  dE 

(E2  -  Az)h 


where  dE  =  de  since  at  T  =  0,  A  is  constant.  Finally  therefore 

IMC«T2D1  (ef)D  (ef)[(eV)2-  A2] 2  2-62 

No  X  L  I  -L 

For  the  case  T  ^  0  the  integral  is  complicated  but  has  been  solved. 

If  we  now  differentiate  this  expression  for  INg  with  respect  to  V 
and  divide  the  result  by  the  derivative  of  1^  with  respect  to  V, 
we  find  that  the  ratio  of  these  dynamic  conductances  is:- 


CTNS 

°NN 


eV 

[(eV)2  -  A2]^ 
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This  expression  is  precisely  the  normalised  density  of  states  for 

a  superconductor  as  can  be  seen  from  the  expression  2—37  derived  earlier. 

We  have  thus  proved  that  by  measuring  the  normalised 

dynamic  conductance  of  an  N  -  S  tunnel  junction  we  are  in  fact 

measuring  the  electronic  quasiparticle  density  of  states  of  the 

superconductor.  Because  of  this,  fundamental  properties  such  as 

T  and  A(0)  are  directly  obtainable, 
c 


7.  Strong  Coupling  Theory 

Whilst  the  weak  coupling  theory  describes  most  superconductors 
extremely  well  and  all  superconductors  to  first  order  there  are 
some  which  deviate  considerably  from  its  predictions.  Mercury  and 
lead  are  two  such  examples.  Experimental  evidence  shows  the 
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condensation  energy  for  lead  to  be  less  than  that  predicted  using 
the  B.C.S.  model  and  outside  the  limits  of  experimental  error. 

The  "universal"  constant  of  equation  2-47  relating  the  gap  to 
transition  temperature  for  lead  was  found  to  be  4.47  instead  of 
3.53. 

Yet  another  striking  difference  is  exhibited  by  these 
superconductors  in  the  density  of  states  as  observed  by  tunnelling 
in  that  there  is  considerable  departure  from  the  smooth  function 
expected  on  the  basis  of  equation  2-37 .  This  was  first  observed  in 
the  tunnelling  experiments  of  Giaver  et.  al.  (1962)  and  Rowell  et.  al. 
(1963)  and  is  illustrated  below  for  the  case  of  lead  which  shows  the 
B.C.S.  prediction  and  the  experimentally  observed  tunnelling  density 
of  states  of  Rowell  et.  al.  (curve  ^2^  )  taken  from  R.D.  Parks. 

t.ia 
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n  i.04 

1.00 

0.86 

0.82 

It  became  clear  in  the  light  of  such  experimental  evidence 
that  the  existing  theory  of  superconductivity  must  be  modified  to 
explain  these  differences  in  behaviour  in  a  class  of  metals  which 
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have  become  known  as  strong-coupling  superconductors.  The  first 
attempt  was  to  try  and  extend  the  weak  coupling  theory  by  removing 
the  assumption  on  the  product  DV.  This  was  done  by  D.J.  Thouless 
(1960)  who  found  that  even  in  the  unphysical  limit  of  DV  00 
the  ratio  2A(0)/kt)T  -*  4.00  and  although  this  was  a  step  in  the 

right  direction  it  could  not  fully  account  for  the  observed 
deviations.  The  next  step  was  taken  by  Swihart  (1962)  who  replaced 
the  simple  B.C.S.  constant  potential  with  a  more  realistic  one. 

This  refinement  met  with  some  success  but  still  could  not  account 
for  the  structure  in  the  density  of  states  curve.  The  correct 
formalism  was  finally  arrived  at  through  the  realisation  that  the 
electron-electron  interaction,  which  according  to  the  B.C.S. 
model  was  instantaneous,  was  in  fact  retarded.  This  can  be  understood 
since  the  phonons  involved  with  the  pairing  of  electrons  move  with 
a  phase  velocity  very  much  less  than  the  Fermi  velocity  of  the 
electrons  and  thus  there  is  a  time  lag  before  the  second  electron 
"feels"  the  existence  (or  is  influenced  by)  the  first.  In  addition 
to  retardation  effects,  lifetime  effects  were  proposed  associated 
with  the  damping  of  quasiparticle  excitations.  This  approach  was 
adopted  by  Migdal  (1958) >  Nambu  and  Eliashberg  (1960)  and  Eliashberg 
(1962)  using  field  theoretic  Greens  functions  (since  potentials 
that  are  retarded  in  time  cannot  be  treated  with  Hamiltonian  dynamics). 


The  results  of  this  analysis  are  presented  below:- 
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where 


t  00 


K*"5 


dva^2(^)f^(v)[(w '  +  oo  +  v  +  i0+)  1±(oo'  -  to  +  v  -  i(T")  1] 


,+s-l 
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Equations  2-64  and  2-65  are  referred  to  as  the  strong-coupling  gap 
equations.  A  third  relationship  is  introduced  in  the  form:- 


A(oo)  =  0  (to  )  /  Z  (to ) 


2-66 


to^  represents  a  cut  off  energy  usually  taken  to  be  -  lOw^,  is 
the  Coulomb  pseudopotential  characterising  the  short  range  repulsive 
force  between  electrons  and  Z(to)  is  a  renormalisation  function  given 


by:- 

e(k,to)  =  e/Z(to) 

where  e  represents  a  single  particle  energy  (Bloch  state). 

We  can  see  that  these  coupled  equations  for  A  are  complex  as  is 
A  itself  in  the  general  case.  Clearly  A  is  also  a  function  of 
energy.  It  may  be  written  in  the  form 

A(co)  =  A  (oo)  +  iA  (to)  2-67 

1  2 

The  a^2 (v)  represents  a  coupling  function  associated  with  the  electron 
phonon  interaction  whilst  the  function  f^(v)  takes  account  of  the 
specific  nature  of  the  phonon  spectrum  of  the  particular  material 
under  investigation.  Thus  we  see  some  marked  differences  in  this 
formalism  from  the  B.C.S.  approach  although  in  the  weak  coupling 
limit  the  quantity  A(oo)  ->■  A(0)  enabling  us  to  identify  A(to)  with 
the  superconducting  energy  gap  still.  To  obtain  A(oo)  we  must  first 
solve  the  integral  equations  2-64  and  2-65  self  consistently  and  then 
apply  equation  2-66.  An  example  of  the  form  of  the  gap  function  is 
shown  below  for  the  case  of  lead  due  to  McMillan  and  Rowell  (1969) 
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Real  ( - )  and  imaginary  (...)  parts  of  the  computed 

energy  gap  function  A(oo)  for  Pb  vs.oo-A(O). 


The  gap  function  assumes  a  constant  value  at  some  given  energy 
which  is  consistent  with  the  minimum  energy  necessary  for  pair 
creation.  The  point  at  which  this  occurs  is  known  as  the  gap 
edge  and  is  the  place  at  (or  below)  which  the  gap  value  is  usually 
measured . 

We  saw  earlier  that  in  the  weak  coupling  approximation  of 
DV  -  10  ^  that  a  universal  constant  of  3.53  resulted  for  the  ratio 
2A(0)/kgTc  and  whilst  this  approximation  is  not  valid  for  strong 
coupling  superconductors  by  definition, it  is  used  instead  as  a  guide 
to  the  strength  of  coupling;  the  higher  its  value  the  stronger  the 
coupling  is  presumed  to  be.  Clearly  since  both  A  and  T^  have  the 
same  functional  dependence  on  the  product  DV,  then  we  would  not 
expect  there  to  be  any  deviation  in  2A(0)/kgTc  as  a  result  of 
increasing  coupling  strength  DV.  Lead  for  example,  however  exhibits 
a  value  of  4.47  for  2A(0)/kgT^.  This  is  thought  to  be  due  to 
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inelastic  phonon  processes  causing  damping  of  quasiparticle  states 

and  although  this  would  reduce  both  A(0)  and  T  ,  the  effect  on 

would  be  greater  since  there  would  be  more  thermal  phonons  present  at 

(as  for  lead  is  high) .  Such  a  calculation  has  been  carried 

out  for  lead  on  this  basis  by  Wada  (1964)  with  a 

prediction  of  2A(0)/kT>T  =  5.2 

B  c 

An  expression  for  2A(0)/kt(T  has  been  derived  for  strong  coupling 

B  C 


superconductors  by  Geilikman  and  Kresin  (1965)  given  by 
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where  co^  is  the  limiting  frequency  of  longitudinal  phonons. 

Another  indicator  as  to  the  strength  of  the  coupling  is  the 
nature  of  the  quantity  DV  =  A  which  is  regarded  as  the  electron-phonon 
coupling  strength.  For  weak  coupling  this  was  assumed  to  be  of  the 
order  10  ^  whilst  for  the  system  under  consideration  here  it  is  of 
the  order  unity.  The  quantity  A  is  given  by 


C  00 


A  =  2 


0 


doxx2  (oj)F  (qj) 


and  1  +  A  =  m* 
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m 

where  m*  is  the  effective  mass  of  the  electron.  A  further  indicator 
is  associated  with  a,  the  coupling  function  appearing  in  the  gap 
equation,  which  is  usually  expressed  as  an  average  given  by:- 


<a2>  = 


/  oo  CO 

a2  (a))F(oo)d w/  F(oo)du) 


2-71 


0;  0J 

Increases  in  the  values  of  A,  <a2>  and  2A  (0)/k  T  are  associated 

’  Be 

with  increases  in  the  coupling  strength.  In  the  experiments  that 
were  conducted  A  and  were  directly  measured  whilst  A  and  <a  > 
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(together  with  other  microscopic  parameters  such  as  y*  the  Coulomb 


pseudopotential)  were  obtained  indirectly  by  a  process  of 
inversion  of  the  strong  coupling  gap  equations.  The  inversion 
program  will  be  discussed  in  more  detail  in  a  later  chapter. 

An  expression  for  the  transition  temperature  of  strong-coupling 
superconductors  has  been  proposed  by  McMillan  (1968)  in  the  form 


T  =  <Ql)>  exp 

C  1.20 


1.04  (1  +  A) 

X  -  y*(l  +  . 62X) 
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If  we  differentiate  this  expression  we  get  (assuming  dy*/dp  =  0):- 


dlnT  =  din  <co>  +1.04 
c 


dp 


dp 


1  +  .38y* 
(A-y*-.62Ay*) 2 


dA 

dp 
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Those  quantities  on  the  right  hand  side  of  this  equation  are 
derivable  from  the  inversion  process  whilst  T^  is  directly  measurable 
and  we  are  thus  able  to  check  the  overall  "quality  of  the  inversion" 
technique  against  experiments. 

Scalapino,  Schrieffer  and  Wilkins  (1963,  1966)  applied  the 
principles  of  strong-coupling  theory  to  the  calculation  of  the  normal 
metal  superconductor  tunnelling  current  and  showed  that  the 


expression  2-63  still  holds  in  the  following  form:- 


a 


NS 

7nn 


Re  eV 


(eV)  -  A2(U)) 


D(E) 

D(£f) 
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Using  a  very  simple  form  for  the  phonon  spectrum  of  lead  which  they 
approximated  by  two  lorrentzians  situated  at  -4.4 meV  and  8.5  meV, 
Schrieffer,  Scalapino  and  Wilkins  (1963)  calculated  the  gap  function  and 
then  the  normalised  conductivity  for  the  NS  case  at  non  zero  temperature. 
These  results  are  shown  in  the  diagram  on  page  30  (fig.  (T)  )  in  which 
we  can  see  that  there  is  excellent  qualitative  agreement  with  the 
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experimental  observations  of  Rowell  et.  al.  (1969). 

8 .  Observation  of  Phonon  Structure  by  Tunnelling 

We  have  now  seen  that  inclusion  of  the  form  of  the  phonon  spectrum 

leads  to  structure  in  the  normalised  density  of  states  curve  obtained 

by  tunnelling  (from  the  results  of  Rowell  et.  al.  for  example). 

That  this  is  plausible  can  be  shown  by  examining  the  form  taken  by 

the  expression  for  the  normalised  conductivity  G  =  in  the 

light  of  the  strong-coupling  theory 

a  =  Re  |  to  1  where  go  =  eV  2-75 

(go2  -  A2  (a>))  *5 

expanding  by  the  binomial  theorem  and  using  equation  2-67  we  get 

a  =  1  +  A  2(aj)  -  A  2  (to)  + .  2-76 

1  _ _ 2 _ 

2G02 

Since  A  (go)  clearly  depends  on  F(go)  (inherent  in  the  strong¬ 
coupling  equations)  then  the  normalised  conductivity  that  we  monitor 
in  electron  tunnelling  experiments  must  also  be  influenced  by  the 
phonon  spectrum  of  the  metal.  Quasiparticles  whose  energies  coincide 
with  the  predominant  phonon  modes  of  a  given  metal  have  a  high 
probability  of  phonon  emission  with  subsequent  decay  to  the  peaked 
region  in  the  density  of  states  (at  the  gap  edge).  We  therefore 
observe  larger  changes  in  the  conductivity  at  these  points. 

Scalapino  and  Anderson  (1963)  have  shown  that  such  resonances  in  F (go) 
give  rise  to  points  of  inflection  in  the  density  of  states.  This  can 

also  be  seen  by  careful  study  of  the  effects  of  A  and  A  on 

12 

using  2-76.  We  are  therefore  provided  with  a  powerful  tool  for 


. 
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resolving  principle  phonon  peaks  of  a  metal  since  differentiation  of 
the  conductivity  curve  will  give  rise  to  maxima  and  minima  at  these 
inflection  points. 

For  the  lead-indium  system  under  investigation  the  principle 
phonon  peaks  are  those  of  the  lead  lattice  but  with  the  addition  of  an 
extra  peak  at  the  high  energy  end  of  the  spectrum.  This  extra  peak 
is  associated  with  inclusion  of  a  small  amount  of  indium  impurity 
in  the  host  lattice.  Justification  for  its  existence  is  the  subject 
of  the  next  section. 


9.  The  Localised  Mode 


To  verify  the  existence  of  a  localised  mode  of  vibration  in  a 
solid  which  contains  a  light  substitutional  impurity  let  us  consider 
the  simple  one  dimensional  case  below :- 
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We  now  consider  the  dynamical  problem  arising  when  one  row  of 

atoms  M  is  replaced  by  the  light  impurity  M*  as  a  longitudinal  wave 

passes  through  the  crystal.  Assuming  only  nearest  neighbour 

interactions  and  further  that  the  force  constant  c  for  these 

interactions  is  constant  (unchanged  by  the  prescence  of  M’)  we 

have  for  the  equations  of  motion  of  the  s  plane  of  atoms 

M’d2u  =  c(u  +  u  -  2u  )  2-77 

_ o  i  -i  o 

dt2 

and  for  the  (s  +  1)  plane 

Md2u  =  c(u  +  u  -  2u  )  2-78 

_ 1  2  0  1 

dt2 


where  the  u  represent  displacements  of  the  planes  of  atoms  from 

their  undisturbed  positions.  Using  a  trial  solution  of  the  form 

,  .  N  s  —  ico t  -I  s  la 
u  =  u  (-1)  e  e  1  1 

s  o 

a  self  consistent  solution  is  obtained  for  2-77  and  2-78  if 
ea  =  2M  -  M' 


M' 


2-79 


hence 


032  =  4c 
M 


M* 


2 MM'  -  M'2 


2-80 


Similar  treatment  of  the  case  when  M  =  M*  (no  impurity)  yields,  for 

„  „  ,  i  (s+p)ka-iojt 

a  travelling  wave  solution  of  the  form  ug  =  u(0)e 

oo2  =  4_c  Sin2ka 
M 

This  has  a  maximum  value  oo2  =  4c/M  thus  substituting  in  2-80  we  get 


2  2 
or,  -i  =a) 
local  max 


M" 


2-81 


2MM'  -  M’2 
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as  the  predicted  frequency  for  the  local  mode  of  vibration. 

For  the  case  at  hand  of  indium  in  lead  we  can  treat  the  indium 
as  such  a  light  impurity.  Indium  (face  centred  tetragonal)  alloys 


with  lead  (face  centred  cubic)  to  concentrations  of  up  to  65  atomic  % 
In  to  give  a  substitutional  f.c.c.  structure.  To  obtain  an 
approximate  value  for  ^  we  put  M  =  207.2,  M*  =  114.8  (the  atomic 
weights  of  Pb  and  In  respectively)  and  take  a)  to  be  the  cut  off 
in  the  Debeye  spectrum  of  lead  and  we  get  that 

E..  _  -  9.83  meV 

local 

This  prediction  is  in  fact  in  excellent  agreement  with  that  observed 
as  will  be  seen.  Whilst  it  is  realised  that  such  agreement  is 
fortuitous  in  the  light  of  the  assumptions  made  it  is  noteworthy 
that  the  more  complex  Greens  function  approach  using  Dysons  equation 
yields  an  expression  identical  to  that  of  equation  2-81. 


CHAPTER  III 


EXPERIMENTAL  APPARATUS 


1 .  The  Cryostat 

A  conventional  double  glass  dewar  system  was  used  to  cool  the 
main  body  of  the  cryostat,  a  cross  section  of  which  is  shown  in 
figure  1.  The  original  design  and  construction  of  the  cryostat  was 
carried  out  by  Franck  and  Keeler  (1968)  for  their  work  on  pure  lead 
and  was  used,  after  extensive  modification,  for  this  work.  The  cryostat 
basically  consisted  of  three  self  contained  systems. 

The  first  of  these  was  a  beryllium-copper  chamber  (which  will 
be  referred  to  hereafter  as  the  bomb)  in  which  the  tunnel  junctions 
to  be  studied  were  actually  mounted.  Commercially  supplied  high 
pressure  capillary  tubing  (Harwood  Eng.  Co.  3M  series)  was  used  to 
connect  the  bomb  to  a  gas  compressor . 

The  second  system  consisted  of  a  small  copper  chamber  (~150  ml) 
attached  to  the  lower  end  of  the  bomb  by  a  large  copper  securing  nut 
and  washer.  Quick  connection  of  this  chamber  to  the  needle  valve  and 
a  roughing  pump  were  made  using  swagelock  couplings  (not  shown  in 
diagram).  By  first  admitting  liquid  helium  into  this  chamber  and 
subsequently  lowering  its  vapour  pressure  by  pumping,  temperatures 
in  the  region  of  IK  -*  2K  could  be  attained.  An  indication  that 
such  temperatures  could  be  realised  can  be  seen  in  figure  2  which 
shows  a  double  gap  in  a  tunnel  characteristic  obtained  from  a  lead- 
insulator-aluminium  junction.  The  presence  of  this  double  gap  is 
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FIGURE  1 
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FIGURE 
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attributed  to  the  onset  of  superconductivity  in  the  aluminium  film 
(T  =  1.2IC)  and  although  it  is  appreciated  that  thin  film  effects  can 
cause  a  rise  in  to  perhaps  2K,  this  is  not  thought  to  be  the  case 
here  since  a  temp  ~1K  was  confirmed  by  the  thermometry.  To  obtain 
the  excellent  thermomechanical  contact  necessary  for  such  low 
temperatures  the  lower  surface  of  the  bomb  was  carefully  lapped  to 
the  upper  surface  of  the  chamber  using  a  series  of  Exolon 
lapping  compounds  and  smearing  the  mated  surfaces  with  a  suspension  of 
copper  powder  in  Apiezon  N  type  vacuum  grease. 

The  third  self  contained  system  took  the  form  of  a  large  copper 

can  (main  vacuum  can)  enclosing  the  other  two  systems  and  connected  to 

a  high  vacuum  oil  diffusion  pump.  Except  when  using  exchange  gas,  the 

—6  -7 

pressure  inside  this  can  was  maintained  in  the  range  10  torr  -*■  10 
torr  and  provided  therefore  the  necessary  thermal  isolation  of  the 
bomb  from  the  cryogenic  fluid  in  which  it  was  situated. 

An  indium  0  ring  seal  was  used  to  retain  the  high  vacuum  at 
the  exit  point  of  the  high  pressure  capillary  from  this  can.  This 
seal  was  located  several  inches  above  the  top  of  the  can  to  provide 
an  extra  degree  of  thermal  isolation  for  the  pressure  tube  and  thus 
reduce  the  possibility  of  blockage  caused  by  cold  spots. 

In  order  to  further  minimise  the  occurrence  of  blockage,  heaters 
were  placed  at  strategic  points  including  a  manganin  resistance  wire 
inside  the  actual  pressure  tube  and  running  its  entire  length. 
Temperature  regulation  was  also  effected  in  part  by  these  heaters. 

Relocation  of  the  needle  valve  from  the  can  top  to  its  present 
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position  permitted  the  use  of  the  very  simple  design  shown  which  was 
not  only  leak  tight  but  since  it  was  always  totally  immersed  in  the 
liquid  helium,  was  not  susceptible  to  icing  up. 

Wires  for  the  heater  circuits  were  admitted  to  the  can  through 
the  high  vacuum  tube  (via  a  room  temperature  epoxy  feedthrough)  and 
thermally  anchored  by  winding  several  inches  of  their  length  on  a 
copper  post  and  securing  them  tightly  with  varnish.  Two  such  thermal 
anchors  were  employed,  one  at  4K  situated  at  the  can  top  and  the 
other  at  ~2K  and  located  at  the  bomb  top.  The  2K  anchor  was  constructed 
in  such  a  way  as  to  be  removable  in  the  event  of  failure  of  the  main 
wiring  (to  be  described) . 

After  passing  through  the  liquid  helium  bath,  the  thermometer 
leads  were  led  into  the  high  vacuum  can  via  an  epoxy  feedthrough  of 
the  type  described  by  Anderson  (1968).  From  here  they  were  then 
connected  to  the  IK  pumping  stage  using  a  jumper  cable  wired  to 
Amphenol  miniplugs  (for  quick  connection)  at  which  point  they  were 
firmly  anchored,  again  using  varnish. 

The  cryostat  was  used  over  a  temperature  range  of  from  IK  40K. 
That  such  elevated  temperatures  were  necessary  can  be  appreciated  from 
the  phase  diagram  (figure  3)  in  which  it  can  be  seen  for  example  that 
in  order  to  form  a  solid  at  3.5  kbar  temperatures  ^35K  are  required. 

2 .  The  Bomb 

Figure  4  shows  a  cross-sectional  viow  of  the  bomb  and  low 
temperature  pumping  stage.  The  bomb  itself  was  machined  from  an  alloy 
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FIGURE  3 
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FIGURE  4 
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containing  approximately  98%  copper  and  2%  beryllium  (with  traces  of 
Co  and  Ni)  which  after  heat  treatment  had  a  tensile  strength  of 
~ 200,000  p.s.i.  The  advantages  of  this  alloy  are  that  it  exhibits 
good  thermal  conductivity  and  has  a  strength  comparable  to  stainless 
steel  without  embrittling  at  low  temperature. 

The  gas  tight  seal  between  the  two  halves  of  the  bomb  was 
achieved  by  compressing  a  specially  designed  stainless  steel  ring 
(lense  ring)  whose  edges  had  been  bevelled  at  37  .  This  particular 
method  provides  a  self  tightening  seal  since  any  increase  in  internal 
pressure  serves  to  further  force  the  ring  onto  the  sealing  edges  of  the 
bomb.  The  core  of  the  bomb  was  machined  to  dia.  and  to  a  depth 

of  ~2^M  to  accommodate  the  sample  housing. 

3 .  Sample  Housing  and  Holder 

In  order  to  maximise  the  useable  space  in  the  bomb’s  experimental 
area  a  special  housing  was  designed  and  constructed  and  is  shown  in 
figure  5  together  with  the  sample  holder . 

The  housing  consisted  of  a  thin  wall  stainless  steel  tube  (whose 
diameter  was  chosen  to  be  a  close  fit  in  the  bomb)  which  had  been 
milled  as  shown.  A  thin  capillary  tube  was  then  swaged  slightly  at 
each  end  and  silver  soldered  to  the  back  of  this  tube  to  allow  passage 
of  the  wires.  An  Amphenol  miniplug  was  next  turned  to  size,  relieved 
so  as  to  accommodate  the  capillary  and  then  epoxied  into  position 
as  illustrated.  Having  attached  a  lid  to  the  top  end,  this  device 
was  then  secured  in  place  in  the  bomb  top.  The  use  of  this  housing 
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allowed  the  passage  of  the  wires  to  the  bottom  end  of  the  plug  where 
they  were  secured  in  such  a  fashion  as  to  be  completely  safeguarded 
against  breakage. 

The  sample  holder  was  constructed  by  first  fabricating  the  base 
from  stainless  steel  and  silver  soldering  to  this  the  beryllium 
copper  clip  arrangement  which  holds  the  glass  substrate  (on  which 
the  tunnel  junctions  are  formed).  Beryllium  copper  was  chosen  for  the 
clips  since  cold  working  of  stainless  steel  leads  to  breakage  in  this 
case.  This  arrangement  was  then  soldered  to  a  compatible  plug.  The 
two  greatest  advantages  of  this  set-up  were  that  it  was  very  fast  to 
mount  a  specimen  whilst  eliminating  completely  the  handling  of  the 
specimen  itself. 

4 .  High  Pressure  Production  and  Measurement 

High  pressure  production  was  achieved  using  a  Harwood  Eng.  Co. 
compressor  specially  adapted  for  use  with  gaseous  helium  four.  From 
an  initial  supply  at  800-1000  p.s.i.  the  gas  was  pumped  in  three 
successive  stages  to  the  3.5  kilobar  range  used  for  these  measurements. 
Isolation  of  the  experimental  chamber  from  the  final  compressor  stage 
was  effected  by  a  needle  valve. 

Three  methods  of  pressure  measurement  were  employed  namely  a 
Heise  Bourdon  gauge  (0-100,000  p.s.i.),  a  manganin  gauge  and  the 
tunnel  junctions  themselves.  The  manganin  gauge  supplied  with  the 
compressor  monitored  the  pressure  dependent  resistance  of  a  manganin 
wire  using  a  Wheatstone  bridge  circuit  and  was  supposedly  accurate 
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to  1%.  This  gauge  and  the  Heise  Bourdon  gauge  were  used  to  calibrate 
the  tunnel  junction  by  recording  its  zero  bias  resistance  as  a  function 
of  pressure.  The  final  frozen  in  pressure  was  always  determined 
using  the  junction  itself  since  the  two  other  methods  would  then  only 
monitor  external  pressure.  This  method  is  to  be  discussed  in  more 
detail  in  a  later  chapter.  Experiment  showed  that  the  temperature 
dependent  resistance  change  of  a  sample  was  <.5%  for  a  AT  of  70K  and 
thus  no  correction  to  pressures  was  applied  in  this  regard. 

5 .  The  High  Pressure  Electrical  Feedthrough 

In  conducting  experiments  at  high  pressures  using  gaseous  helium 
as  the  pressure  transmitter  one  is  faced  with  the  problem  of 
introducing  electrical  wires  into  the  experimental  chamber  without 
loss  of  pressure.  To  overcome  this,  a  seal  was  constructed  based 
on  the  frozen  oil  method  employed  by  Dugdale  and  Hulbert  (1957)  and 
is  illustrated  in  figure  6.  An  oil-epoxy  feedthrough  which  incorporates 
this  principle  has  been  described  by  Goree  et.  al.  (1965),  in  which 
heavy  emphasis  is  placed  on  the  role  of  the  epoxy  which  separates  the 
gas-oil  interface.  This  particular  seal  was  found  to  be  extremely 
difficult  to  construct  and  its  failure  rate  in  service  found  to  be 
high,  both  of  these  drawbacks  being  attributed  to  the  presence  of 
the  epoxy.  Furthermore  a  close  examination  of  the  defunct  feedthroughs 
revealed  that  the  epoxy  had  in  fact  not  penetrated  the  pressure  tube 
to  any  appreciable  extent.  This  latter  problem  was  due  in  part  to  the 
fact  that  in  the  experiments  described  herein,  6-9  wires  were  necessary 
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FIGURE  6 
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thus  providing  a  major  constriction  in  the  pressure  tubing. 

To  this  end,  a  feedthrough  was  devised  for  use  with  large 
numbers  of  wires  which  did  not  employ  the  use  of  epoxy  filled  tubes. 
This  feedthrough  is  in  comparison  very  easy  to  construct  and  re— 
tillable  as  well  as  incorporating  a  plug  designed  to  protect  wires  from 
accidental  damage. 

Subject  to  a  specifically  ordered  filling  procedure  which 
prevents  contamination  of  the  experimental  chamber ,  the  feedthrough 
illustrated  not  only  survived  the  test  of  time  but  vented  itself 
of  oil  when  the  liquid  nitrogen  level  fell;  again  without  contamination 
of  the  experimental  chamber. 

In  order  to  construct  the  particular  seal  shown,  approximately 
twelve  feet  of  3/16M  O.D. ,  .0025"  bore  stainless  steel  pressure 
tubing  (Harwood  3M  series)  was  first  coned  and  threaded  (for  use 
with  the  standard  fitting  shown)  and  cut  to  the  desired  lengths  for 
the  tubes  T1 ,  T2 ,  and  T3.  These  tubes  were  then  laid  out  on  a  long 
clean  surface  together  with  the  junction  blocks  and  fittings  in  a 
development  which  would  ultimately  lead  to  the  final  form  shown. 

Having  passed  the  required  number  of  wires  through  the  system  and 
allowing  an  excess  at  each  end,  the  fittings  were  carefully  tightened 
and  the  feedthrough  bent  to  shape. 

The  principle  involved  in  retaining  the  gas  pressure  with  such 
an  arrangement  would  be  to  fill  the  tube  T3  with  oil  by  admitting  it 
through  0  in  block  C  and  then  immersing  T3  in  liquid  nitrogen.  In 
order  to  prevent  the  preferential  leakage  of  oil  at  the  top  of  block 


53 


C  during  filling  it  was  necessary  to  block  T4.  To  this  end  a  re¬ 
usable  device  (upper  seal  in  illustration)  was  constructed  which 
not  only  blocked  T4  but  accommodated  a  plug  so  that  the  wiring 
could  be  accessed  without  exposure  to  possible  breakage.  Since  T4 
was  never  called  upon  to  withstand  more  than  two  or  three  thousand 
p.s.i.  oil  pressure,  a  slightly  larger  —  but  still  standard  — 
pressure  tube  was  used  having  an  o.d.  of  1/4"  and  a  bore  of  1/16 
(Harwood  4L  series) . 

To  assemble  the  upper  seal,  the  main  body  of  machined  brass  (U3)  and  the 
the  stainless  steel  cup  (U5)  were  first  slid  onto  the  tube  followed 
by  the  standard  threaded  collar  (U6).  The  stainless  steel  cup  was 
held  in  position  and  filled  with  Stycast  2850  G.T.  epoxy  which  was 
allowed  to  cure.  Some  of  this  epoxy  is  also  expected  to  enter  T4  and, 
although  favourable,  is  not  essential  to  the  operation  of  the  seal. 

The  protruding  wires  were  then  attached  to  an  Amphenol  miniplug  (Ul) 
housed  in  the  machined  brass  cap  (U2).  Finally  the  main  body  was 
raised  into  position  and  locked  onto  the  tube  and  the  cap  secured. 

To  fill  the  seal,  very  low  pressure  helium  gas  was  pumped  into 
block  A  through  G  and  allowed  to  leak  slowly  from  the  vent  V. 

Simultaneously,  oil  was  pumped  at  low  pressure  into  block  C  through 
0  and  thus  through  tube  T3  to  block  B  and  also  allowed  to  leak  from 
V.  It  was  felt  that  the  steady  stream  of  gas  in  T2  would  prevent 
the  slight  possibility  of  oil  leaking  into  this  tube.  Having 
established  a  continuous  bubble  free  flow  of  oil  from  the  vent,  the 
tube  T3  was  immersed  in  liquid  nitrogen.  When  the  oil  ceased  to 
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flow,  the  oil  pressure  was  reduced  to  zero  and  the  vent  valve  V 
tightened  whereupon  the  system  pressure  was  observed  to  rise  as 
expected.  The  bomb  remained  in  position  throughout  the  filling 
process.  The  system  described  is  now  ready  for  use  and  will  remain 
so  if  the  nitrogen  level  is  maintained. 

When  not  in  active  use  (or  when  deliberately  venting  the  oil 
seal)  the  bomb  was  kept  in  place  containing  a  few  hundred  p.s.i. 
of  gas  pressure  while  the  oil  supply  line  0  remained  vented.  This 
practice  ensured  that  if  the  nitrogen  trap  was  left  untended  the 
melted  oil  would  be  ejected  from  0  and  eliminate  any  possibility 
of  contamination  of  Tl  via  T2.  Should  this  occur  it  would  be  necessary 
to  recharge  the  seal  by  repeating  the  above  procedure. 

6 .  Vacuum  Deposition 

The  tunnel  junctions  were  formed  by  evaporation  in  a  customised 
Norton  Research  Company  3115  model  high  vacuum  system.  Masks  were 
fashioned  from  stainless  steel  to  give  the  required  film  geometry 
and  an  enlarged  contact  area  and  located  in  a  carousel  which  allowed 
up  to  six  different  source-mask  combinations  to  be  used.  The  carousel 
was  enclosed  in  a  large  glass  tube,  with  quartz  glass  divisions  to 
prevent  mutual  contamination  of  sources  and  masks. 

An  arrangement  compatible  with  the  sample  holder  described 
previously  was  incorporated  in  the  carousel  design  in  such  a  way  as 
to  allow  any  desired  mask  to  be  raised  into  contact  with  the  substrate. 
Since  the  back  of  the  sample  holder  had  been  machined  away,  it  was 
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possible  (with  the  aid  of  a  prism)  to  view  the  film  during  formation. 

7 .  Temperature  Measurement 

Temperatures  were  measured  with  the  aid  of  a  commercially  supplied 
calibrated  germanium  resistance  thermometer  (Cryocal),  the  calibration 
range  being  from  1.5K  to  40K. 

The  1.5K  to  2. OK  readings  were  based  on  the  Helium  Vapor  Pressure 
Scale  of  1958  with  an  accuracy  of  0.003K  while  the  2.25K  to  13K 
readings  are  traceable  to  the  National  Bureau  of  Standards  Provisional 
Scale  of  1965  (2-20°K)  through  secondary  standards  maintained  by 
CryoCal,  Inc.  with  an  accuracy  of  0.005K  from  2K  to  5K  and  0.01K  from 
5K  to  13K. 

The  14K  to  40K  readings  are  traceable  to  the  International 
Practical  Temperature  Scale  of  1968,  again  through  secondary  standards 
maintained  by  CryoCal,  Inc.  with  an  accuracy  of  0.01K  from  14K  to 
20K  and  0.04K  from  22K  to  40K. 

Resistance  values  were  obtained  to  better  than  .14  using 

potentiometric  techniques. 

Standard  curve  fitting  procedures  were  employed  for  interpolation 
(and  extrapolation  where  necessary)  with  excellent  reliability. 


CHAPTER  4 


ELECTRONIC  DETECTION  TECHNIQUES 

1.  The  Principle  of  the  Measurement 

It  was  shown  earlier  (Chapter  2)  that  the  quasiparticle  density 
of  states  in  a  superconductor  was  in  one  to  one  correspondence  with 
the  normalised  dynamic  conductance  C^g/a^  as  measured  in  a  normal- 
superconductor  (NS)  tunnel  junction.  Measurements  of  this  quantity  were 
made  using  standard  harmonic  detection  techniques  by  applying  a  small 
constant  a.c.  signal  across  the  tunnel  junction  and  monitoring  its 
a.c.  response.  Considering  the  idealistically  simple  circuit  shown 
below 
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where  a  constant  a.c.  signal  6Cosiot  (current  or  voltage)  is  applied 
across  the  sample  R^  together  with  a  d.c.  bias  current  Iq.  The 
signal  developed  across  the  sample  resistance  R^  therefore  may  be 
expressed  using  a  Taylor  series  expansion  as:- 


56 


57 


V(I)  =  V(I  +  SCoscot)  =  V (I  )  + 
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where  6Coscot  here  represents  a  constant  a.c.  current  supply. 

If  we  now  let  6  Cosoot  represent  instead  a  constant  voltage  supply 


we  obtain  in  a  similar  manner :- 


I(V)  =  I(V  +  6 Cosoot)  =  I(V  )  + 
O  o 


dl 
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where  in  both  equations  4-1  and  4-2  use  has  been  made  of  a  trigonometric 
identity.  Thus  we  can  see  that  if  we  can  measure  the  components  of  the 
signal  across  R  that  are  proportional  to  0)  and  2to  we  are  in  fact 
measuring  quantities  directly  proportional  to  the  dynamic  resistance 
(or  conductance)  and  its  derivative.  The  component  proportional  to  to 
is  used  to  yield  a  /a  whilst  the  component  proportional  to  2to 

IN  O  IN  IN 

is  used  to  resolve  structure  in  the  density  of  states. 


2 .  Detection  of  the  Harmonics 

In  order  to  measure  the  amplitudes  of  the  components  of  to  and 
2 co,  phase  sensitive  detection  techniques  were  employed  using  the 

circuit  logic  shown  in  figure  7. 

An  a.c.  signal  is  applied  across  the  tunnel  junction  which  is 
incorporated  as  one  arm  of  a  bridge  circuit.  By  so  doing  we  are 
able  to  balance  out  the  large  component  of  fundamental  and  observe 
only  the  small  non  linear  changes  associated  with  dV/dl  etc.  The 
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FIGURE 
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response  of  the  tunnel  junction  is  then  fed  to  the  input  of  a  phase 
sensitive  lock-in  amplifier  (P.A.R.  HR8)  whose  phase  is  set  initially 
to  zero  and  whose  band  pass  is  centred  on  the  source  frequency. 
Simultaneously  a  reference  signal  from  the  same  source  is  fed  to  the 
lock— in  amplifier  and  in  this  way  we  are  able  to  observe  the 
component  proportional  to  oj.  The  source  of  a.c.  supply  should  be  as 
stable  in  amplitude  as  possible  and,  less  critically,  stable  in 
frequency.  Variations  in  oscillator  amplitude  are  most  serious 
since  not  only  do  they  manifest  themselves  directly  at  the  lock  in 
output}  but  sweep  times  associated  with  each  curve  traced  might  range 
up  to  an  hour  in  duration.  Slight  drift  in  frequency  is  not  so  serious 
(depending  on  the  band  pass)  since  the  reference  supply  is  derived 
from  the  source.  The  bridge  circuit  is  balanced  by  minimising  the 
output  from  the  lock-in  amplifier  (for  a  specific  d.c.  bias  condition) 
both  resistively  and  capacitively ,  in  the  latter  case  the  phase  at 
the  lock-in  is  shifted  by  90°.  Having  done  this  (and  returned  the 
phase  of  the  detector  to  zero)  unbalancing  of  the  capacitive  circuit 
should  cause  no  appreciable  effect  on  the  output  of  the  lock  in. 

Finally  while  locked  onto  the  resistive  component  of  the  sample  the 
d.c.  bias  is  changed  slightly  causing  an  unbalanced  bridge  condition 
and  corresponding  output  signal  from  the  lock-in.  Fine  adjustments 
of  the  phase  and  frequency  are  then  made  to  maximise  this  output. 

The  system  is  thus  tuned  for  maximum  response  to  the  first 
harmonic  content.  In  order  to  detect  the  second  harmonic  component 
it  is  only  necessary  to  change  the  oscillator  frequency  to  w/2.  while 
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leaving  the  lock-in  tuned  to  u). 

The  output  from  the  lock-in  is  then  fed  to  a  recorder  where 
it  is  plotted  as  a  function  of  applied  d.c.  bias  using  an  X  -  Y 
recorder.  The  desired  amplification  required  is  achieved  by  the 
combined  effect  of  the  variable  gains  of  lock-in  and  pen  recorder. 

3.  Circuit  Design  Considerations 

Since  the  tunnelling  experiments  carried  out  were  in  the  Fermi 
surface  region  it  is  essential  that  the  measuring  technique  itself 
should  not  interfere  with  the  Fermi  surface  shape.  The  temperature 
dependence  of  the  Fermi  surface  shape  is  often  alluded  to  as  thermal 
smearing,  an  effect  which  becomes  less  pronounced  as  the  temperature 
decreases.  At  IK  for  example  the  thermal  smearing  is  -  86vieV  r.m.s. 
To  this  end  it  is  essential  that  the  a.c.  level  (modulation) 
applied  across  the  tunnel  junction  by  the  circuit  be  as  small  as 
possible  and  less  than  the  thermal  smearing  at  the  temperature  of 
the  measurements. 

In  the  expansions  given  by  equations  4-1  and  4-2,  the  component 
of  the  first  harmonic  is  already  small  compared  to  that  of  the 
fundamental  but  with  the  further  restriction  of  low  applied  signal 
levels  and  the  fact  that  the  tunnel  junction  response  is  only  weakly 
non  linear,  it  is  necessary  to  incorporate  the  tunnel  junction  in  a 
bridge  circuit  in  which  the  large  fundamental  component  can  be 
balanced  out. 

Two  such  bridge  circuits  were  employed  for  this  purpose  one 
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being  a  constant  current  device  of  the  type  used  by  Adler  and  Jackson 
(1966)  and  the  other  a  constant  voltage  device  due  to  Rogers  (1970). 
These  circuits  are  shown  in  detail  in  Figures  8  and  9  respectively, 
the  former  showing  the  support  systems  which  were  used  for  both. 

The  constant  voltage  bridge  was  used  for  measurements  in  the 
immediate  energy  gap  region  where  the  sample  resistance  approaches  °°. 
This  avoids  the  problem  of  over  modulation  associated  with  the 
constant  current  bridge  and  consequent  smearing  of  data.  Both  bridge 
circuits  allowed  four  terminal  measurements  to  be  made  and  included 
large  resistors  to  swamp  the  normal  state  resistance  of  the  film 
electrodes.  The  detailed  circuit  for  d.c.  biasing  (sweep)  is  shown 
in  figure  10. 

4.  The  Constant  Current  Bridge 

The  diagram  below  shows  a  simplified  form  of  the  resistance  bridge 
circuit  used  in  which  R^  represents  the  sample  and  R^  the  bridge 
balancing  components. 


d.c.  bias  (sweep) 


CONSTANT  CURRENT  BRIDGE  (  and  support  systems  ) 
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If  the  current  in  each  arm  of  the  bridge 
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=  at  balance.  Thus:- 
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We  see  therefore  that  detecting'  "^n  ^ac*"  measuri-n§  the 

required  quantity  dV/dl  =  a 

Furthermore  it  can  be  seen  that  the  bridge  response  is  linear  as 
would  be  expected  in  a  passive  device.  This  bridge  is  claimed 
by  its  authors  to  be  able  to  measure  resistance  changes  of  a  few 
parts  in  10~*. 


Calibration  of  the  characteristics  obtained  can  be  achieved 


by  replacing  1^,  by  a  standard  variable  resistance. 


5.  The  Constant  Voltage  Bridge 

This  circuit  is  of  the  active  variety  and  actually  measures 
conductance  as  opposed  to  resistance.  A  high  gain  operational 
amplifier  is  used  to  provide  the  necessary  inversion  to  do  this. 

A  fraction  of  the  output  voltage,  a  say,  is  fed  back  so  as  to  maintain 
a  zero  signal  condition  at  the  input  subject  to  the  condition  that 
aA»l  where  A  is  the  gain. 

Conductance  measurements  of  a  few  parts  in  10  are  obtainable  to 
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an  accuracy  of  1%  and  the  response  of  the  circuit  is  linear  to  .1%. 
Since  the  circuit  contains  an  active  component  a  proof  of  its  linear 
response  is  given. 

The  constant  voltage  bridge  may  be  represented  in  simplified 
form  as  shown  below. 


By  definition 

eout  =  AelN  =  A(e+  '  e_)  =  A(IRD  -  IRT  '  5IV 

According  to  the  dotted  path  to  ground  we  can  see  that 

e  =  6lRp  +  (I  +  6l)Rp  =  <5l  (Rp  +  Rp)  +  IRp 


equating  4-3  and  4-4  gives  after  simplification 
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5l  =  I 


l^Ol-A)  +  ARd 


4-5 


Rp  +  Rj,  +  AR^ 

v. 

Extracting  a  factor  3  =  R,p/(Rp  +  ^  we  can  wri-te  4-5  as:- 
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If  the  feedback  resistor  R„  had  a  value  of  10  ohm  with  the 

r 

amplifier  gain  being  -  10^  then  for  tunnel  junction  resistances 
R^,  >,  10  the  prefactor  in  equation  4-6  is  approximately  unity.  Thus 
neglecting  the  small  quantity  1/A  and  writing  the  tunnel  conductance 


1/R^,  as  Gp  we  have 

51  =  itRjjGp,  -  1] 
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substituting  into  equation  4—4  and  simplifying  yields:- 

eout  ■  (IRdVGT  +  I(RD  -  V  4‘8 

Clearly  therefore  the  output  of  this  bridge  is  also  linear  in  response 
but  this  time  in  conductance.  To  calibrate  the  curves  obtained,  two 
points  only  are  required.  The  most  convenient  points  being  those 
obtained  in  an  open  circuit  condition  (zero  conductance)  and  at  the 
balance  condition  where  G^=  1/R^. 

The  condition  imposed  upon  allowable  sample  resistances  by  Rp 
and  A  were  well  satisfied  in  the  experiments  performed. 


CHAPTER  5 


EXPERIMENTAL  METHOD 


1 .  Sample  Fabrication 

Prior  to  the  fabrication  of  the  tunnel  junction  itself,  alloys 
containing  various  known  atomic  concentrations  of  indium  in  lead 
were  produced  in  a  cleaned  environment  from  extremely  high  purity 
constituents  (69  purity)  which  were  supplied  commercially  (Cominco 
Co.)*  The  tunnel  junctions  formed  from  these  alloys  were  of  the 
traditional  metal-insulator-metal  type  prepared  by  vacuum  deposition. 
Aluminium  was  chosen  for  the  second  electrode  for  two  main  reasons. 
First,  it  does  not  itself  exhibit  any  appreciable  structure  which 
might  cloud  observation  of  the  alloy  in  the  energy  range  of  interest 
and  second,  it  oxidises  readily  thereby  providing  a  very  convenient 
method  of  producing  the  insulating  barrier. 

A  1"  by  1/4"  section  was  cut  from  a  standard  glass  microscope 
slide  and  after  being  cleaned  and  fire  polished  was  mounted  in  the 
sample  holder  previously  described.  This  provided  the  substrate  for 
the  tunnel  junction  and  from  this  point  onward  it  was  never  necessary 
to  touch  the  slide  again.  The  sample  holder  was  then  constrained  in 
a  jig  while  indium  metal  contacts  were  applied  to  the  surface  of  the 
glass  in  places  predetermined  by  the  evaporator  geometry.  Six  such 
contacts  were  made  and  wired  to  the  amphenol  plug  in  order  to  allow 
four  terminal  measurements  to  be  made  on  two  junctions  simultaneously. 

The  sample  holder  was  placed  in  the  evaporator  which  was 
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pumped  to  a  pressure  of  better  than  1  x  10  torr  before  depositing 
the  aluminium  electrode  from  a  multistranded  tungsten  filament. 

The  sample  was  then  removed  from  the  evaporator  and  placed  in  a 
Fisher  Iso temp  oven  in  order  to  form  an  oxide  layer  on  the 
aluminium.  Oxidation  temperatures  used  ranged  from  100°C  to 
130°C  and  were  applied  for  durations  of  from  1  minute  to  5  minutes. 

This  resulted  in  final  tunnel  resistances  of  from  50  ohm  to  300  ohm. 
Several  attempts  were  made  to  control  the  oxidation  process  within 
the  evaporator,  one  of  these  consisting  of  a  glow  discharge  in  low 
pressure  oxygen  and  while  the  glow  discharge  method  was  very 
reproducible  the  oxide  layers  so  produced  did  not  in  general  survive  the 
application  of  pressure.  Having  returned  the  sample  to  the  evaporator, 
the  alloy  electrodes  were  deposited  using  a  molybdenum  boat.  Both 
the  alloy  and  aluminium  electrodes  had  an  approximate  width  of  1/16" 
with  the  latter  running  the  length  of  the  substrate. 

2 .  Mounting  Procedure 

After  filling  the  oil  seal,  a  slow  steady  stream  of  helium  gas 
was  allowed  to  pass  through  the  compressor  to  the  experimental  area. 
Simultaneously  the  electronic  detection  system  was  readied  for  an 
immediate  check  on  tunnel  resistance.  The  sample  was  then  removed 
from  the  evaporator  and,  subject  to  a  satisfactory  visual  inspection 
of  the  films,  was  plugged  into  the  specially  designed  housing 
whereupon  it  was  possible  to  check  instantly  the  approximate  tunnel 
resistance.  This  was  achieved  by  presetting  the  modulation  level  for 
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a  nominal  100  ohm  sample  and  noting  the  change  upon  sample  contact. 

The  bomb  was  then  screwed  loosely  into  place  while  being  flushed 
with  helium  during  which  time  a  more  accurate  check  of  the  tunnel 
resistance  was  made. 

After  tightening  the  bomb  into  place  the  helium  pumping 
chamber  was  secured  (as  shown  in  figure  4)  and  thermometer 
connections  established.  The  vacuum  can  was  then  soldered  into 
place  using  Woods  metal  and  evacuated  to  a  pressure  of  approximately 
lO-'*  torr  at  which  point  the  system  was  cooled  to  liquid  nitrogen 
temperature.  Pumping  of  the  system  was  then  continued  to  the 
lowest  attainable  pressure  in  order  to  provide  good  thermal 
isolation  between  the  cryogenic  fluid  and  the  bomb. 

3.  Validity  Checks 

Before  collecting  all  of  the  data  required  from  a  sample, 
certain  checks  were  made  both  on  the  sample  itself  and  the  system. 

At  liquid  nitrogen  temperature  a  tunnel  characteristic  was  recorded 
for  E  =  +  30  meV  to  E  ~  -  30  meV  in  order  to  check  for  major 
abnormalities  and  noise  levels.  Pressure  approaching  that  used  for 
actual  measurements  was  then  applied  while  monitoring  the  zero  bias 
resistance  of  the  sample  to  establish  reproducibility  in  this  respect. 
In  addition  this  provided  a  check  on  the  pressure  system  itself 
including  the  feedthrough.  Between  the  coo]  down  from  liquid  nitrogen 
to  liquid  helium  temperatures,  the  needle  valve  operation  was  checked. 
The  transition  temperature  was  also  recorded  during  this  process. 
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The  sample  was  then  further  cooled  to  the  ultimate  limit  attainable 
and  examined  for  the  existence  of  the  localised  mode.  This  was 
most  easily  resolved  by  observing  the  second  derivative  curve 
d2V/dI2.  At  this  temperature  the  resistance  in  the  gap  region  at 
zero  bias  was  checked  in  comparison  to  its  normal  state  value. 

Since  by  requiring  that  this  ratio  be  <10  3  at  ~1K  we  can  safely  neglect 
any  non- tunnel ling  contribution  to  the  current. 

4.  An  Experimental  Run 

Having  established  the  validity  of  a  given  sample  and  that 
support  systems  were  operating  efficiently  the  various  low  temperature 
zero  pressure  characteristics  of  the  sample  were  recorded.  The 
cryostat  was  then  left  until  the  level  of  the  He^  bath  was  just 
below  that  of  the  vacuum  can.  It  was  usually  found  that  this  practice 
resulted  in  a  temperature  slightly  in  excess  of  that  predicted  by 
the  Simon  equation  (1937  and  1953)  for  working  pressures  of 
approximately  50,000  p.s.i.g.  During  the  actual  pressur isation 
process,  heaters  around  and  inside  the  pressure  tubing  were 
used  to  insure  against  capillary  blockage.  Monitoring  of  the  zero 
bias  tunnel  resistance  as  a  function  of  pressure  was  also  carried 
out  at  this  time  for  the  prime  purpose  of  determining  the  final 
frozen  in  pressure  although  it  was  also  possible  to  detect  the 
occurrence  of  capillary  blockage  since  when  this  occurred  the 
external  pressure  rose  while  the  sample  resistance  remained  constant. 

Freezing  at  constant  pressure  was  then  achieved  by  opening 
the  needle  valve  while  simultaneously  maintaining  the  external  pressure 
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at  a  constant  value.  In  view  of  the  location  of  the  IK  cooling 
stage  it  was  expected  that  the  solid  would  form  from  the  bottom 
of  the  bomb  upwards. 

On  the  occasions  when  both  junctions  on  the  same  substrate 
were  useable  it  was  possible  to  determine  the  pressure  at  two 
points  in  the  solid  separated  by  -  V'.  The  pressure  variation  over 
this  distance  was  found  to  be  of  the  order  of  1%  and  greatly  supports 
the  belief  that  the  technique  used  did  in  fact  produce  pressures 
which  were  hydrostatic.  With  the  sample  now  suspended  in  solid 
helium  at  constant  pressure,  the  liquid  helium  was  restored  to  a 
level  several  inches  above  the  indium  0  ring  seal  and  the  cryostat 
allowed  to  cool  to  the  region  of  IK,  at  which  point  the  finite 
pressure  characteristics  were  recorded. 

The  penultimate  stage  of  each  experimental  run  provided  the 
challenge  of  releasing  the  pressure  without  destroying  the  sample 
and  the  wiring.  It  was  naturally  not  possible  to  attain  the 
temperature  necessary  to  do  this  while  the  level  of  the  helium  was 
high  and  so  heaters  were  employed  to  boil  off  the  helium  to  a  level 
below  the  can.  Heat  was  then  applied  to  the  bomb  and  the  high 
pressure  capillary.  When  the  temperature  reached  a  value  just  below 
that  at  which  freezing  in  occurred,  the  external  pressure  was 
raised  to  just  above  the  frozen  in  pressure  and  heating  continued. 
When  the  fluid  phase  was  realised  the  pressure  indicated  by  the 
Heise  Bourdon  gauge  showed  a  sudden  change  in  value  corresponding 
to  the  equalisation  of  pressures.  This  change  was  also  noted  at 
the  sample  itself  via  the  zero  bias  resistance.  Early  attempts  to 
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release  the  pressure  at  temperatures  slightly  above  the  melting 
temperature  resulted  in  destruction  of  both  sample  and  wiring. 

This  was  thought  to  be  due  to  the  retention  of  a  plug  of  solid 
helium  high  in  the  capillary  tube  which  itself  is  evidence  that 
the  cryostat  temperature  is  non  uniform  during  this  process. 

An  alternative  method  of  de-pressurisation  would  be  to  wait  until 

the  helium  level  fell  naturally  to  the  desired  level  but  since 

the  cryostat  must  be  full  duing  the  actual  measurements,  this 

would  always  be  an  extremely  lengthy  wait.  In  contrast,  the 

method  used  involved  the  use  of  elevated  temperatures,  a  disadvantage 

which  was  tolerated  subject  to  the  careful  monitoring  of  the  process. 

After  returning  the  system  pressure  to  zero  again  a  further 
check  on  the  reproducibility  was  made  both  on  the  resistance  value  and 
the  observed  structure  location. 

5 .  Determination  of  Sample  Composition 

Despite  the  fact  that  the  sample  was  formed  from  an  alloy 
whose  composition  was  known,  the  final  composition  of  the  evaporated 
electrode  remained  to  be  determined  since  the  vapour  pressures  of  the 
elements  differ  widely  (see  for  example  Nesmayanov  (1963)).  To  do 
this,  residual  resistance  ratio  techniques  were  employed  (see  for 

example  Rosenberg)  in  which  the  ratio  R293/R4  f°r  the  samples  was 
compared  to  a  calibration  curve  of  the  same  quantity  that  had  been 
obtained  by  measuring  a  series  of  alloys  of  known  composition. 

Several  pieces  of  each  known  alloy  were  taken  from  various  places 
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on  the  ingot  and  rolled  out  to  -  .001"  in  thickness.  The  same 
sample  holder  that  was  used  for  the  tunnelling  measurements  (figure  5) 
was  used  to  mount  these  speciments  but  with  a  slight  variation.  A 
plain  glass  slide  was  first  placed  in  the  holder  and  four  indium 
contacts  were  soldered  across  its  width.  Wires  were  then  attached 
to  these  contacts  at  the  edges  and  also  to  the  plug.  Pieces  of  the 
rolled  ingot  were  then  cut  to  slightly  less  than  the  size  of  this 
glass  slide  and  placed  on  top  in  such  a  way  as  to  cover  the  raised 
indium  contacts.  A  second  slightly  smaller  glass  slide  was  then 
placed  on  top  (thereby  sandwiching  the  sample  onto  the  contacts)  and 
held  in  place  with  masking  tape. 

The  sample  holder  was  then  placed  in  a  dipstick  arrangement 
wired  for  a  four  terminal  resistance  measurement.  Resistance 
measurements  were  then  made  using  a  Cambridge  Instruments  high 
resolution  potentiometer  both  at  room  temperature  and  at  liquid 
helium  temperature.  The  latter  measurement  was  made  by  placing  the 
sample  at  the  centre  of  a  Westinghouse  21  Kilogauss  superconducting 
magnet  and  applying  the  necessary  field  to  induce  the  normal  state. 
All  measurements  were  carried  out  for  currents  in  both  directions 
and  averaged  over  all  samples  from  an  ingot  of  a  particular 
concentration.  It  was  in  fact  found  in  this  way  that  the  original 
alloys  were  homogeneous.  The  calibration  curve  obtained  by  this 
procedure  is  shown  in  figure  11. 

A  determination  of  the  ratio  ^293^4  ^or  Pure  ^ea<^  from 
which  the  alloys  were  made  yielded  a  value  of  -  21,000  thereby 


FIGURE  11 

R(293)/  R(4)  Vs  Indium  Concentration 
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attesting  to  its  high  degree  of  purity.  Before  the  ratio  ^293^4 
for  the  actual  films  that  were  tunnelled  through  could  be  determined 
it  was  necessary  to  apply  two  more  contacts  to  the  film  to  enable 
a  four  terminal  resistance  measurement.  This  was  done  as  shown  in 
the  diagram  below  which  shows  an  enlarged  view  of  one  tunnel 
junction. 


Aluminium 


Pb  In, 
x  1-x 

Alloy 


Glass  Substrate 


In  Solder  Contacts 
Ag  Paste  Contacts 


Indium  was  first  deposited  at  two  places  as  close  as  possible  to  the 
existing  contacts  and  wires  attached.  Bridges  between  these  points  and  the 
film  itself  were  then  formed  using  silver  paste.  Several  such  coats 
of  silver  paste  resulted  in  very  low  resistances  from  A  -  B  of  the 
order  of  a  fraction  of  an  ohm.  Once  again  the  same  sample  holder 
was  utilised.  Indium  metal  was  chosen  throughout  for  contacting 
purposes  due  to  its  ability  to  "wet"  glass,  however,  for  this 
particular  application  it  was  essential  not  to  attempt  to  solder 
to  the  film  itself  with  indium  since  not  only  would  the  film  probably 
not  survive  the  process  but  more  importantly  if  it  did,  it  would  not 
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be  possible  to  say  with  any  certainty  that  indium  from  this  source 
had  not  itself  diffused  into  the  film.  Contacts  using  silver  paste 
alone  to  secure  the  wires  were  tried  but  found  to  be  unreliable 
at  low  temperature. 

It  was  very  reassuring  to  notice  during  these  experiments  that 
the  potential  difference  across  the  alloy  films  vanished  (as  would 
be  expected)  as  the  superconducting  state  was  realised  regardless 
of  the  nature  of  the  contacts. 


CHAPTER  VI 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 


1.  Introduction 


It  has  been  observed  by  several  workers  (see  Adler  (1967)  for 
example  and  Sood  (1971))  that  the  major  properties  of  interest  do  not 
differ  significantly  over  samples  whose  impurity  concentrations  vary 
by  ~2  at  %  which,  coupled  with  the  inability  of  a  large  proportion 
of  the  samples  to  withstand  high  pressure,  led  this  author  to  limit 
investigations  to  the  concentrations  Pb^In^  and  P^gg^11^’  Such  a 
choice  offers  the  advantage  of  being  able  to  group  together  data 
from  samples  whose  impurity  concentrations  vary  by  the  above  mentioned 
tolerance  (2%)  while  still  being  indicative  of  any  trends  in  the 
data  which  are  concentration  dependent  (as  a  consequence  of  the 
~25  at  %  difference).  Except  where  specifically  mentioned  therefore, 
results  quoted  henceforth  indicate  averages  for  the  particular 
compositions  over  those  samples  whose  concentrations  (as  determined 
by  the  residual  resistance  ratio  methods  described  earlier)  satisfied 
the  above  condition. 

The  capacitance  of  each  sample  was  determined  in  order  to 
estimate  the  thickness  of  the  oxide  layer  (which  constituted  the 


barrier).  These  values  ranged  from  0.038yF  to  0.051yF  over  nine 
samples  which  included  both  the  Pbggln^  an<3  Pb^In^  samPles  (since 
the  oxide  layer  is  dependent  only  on  the  oxidation  conditions  of 
the  aluminium  electrode  which  is  common  to  both  alloys).  Assuming 
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the  tunnel  junctions  to  be  simple  parallel  plate  capacitors  of  area 

O 

2.25  x  10  m  and  obeying  the  relationship  C  =  £A/d,  the  barrier 

o 

thickness  is  estimated  to  be  ~50A  where  the  dielectric  constant 
of  the  A^O^  has  been  taken  as  9.34  from  Gray  (1972). 

By  measuring  the  resistance  of  the  alloy  electrodes  (which  were 
of  the  order  1.5  ohm  -  2  ohm)  and  assuming  the  value  of  21yohm  cm 
for  the  resistivity  of  lead,  as  taken  from  Meaden  (1965),  the 

0 

thickness  of  the  alloy  electrodes  is  estimated  to  be  from  4400A  - 
0 

5800A. 

2 .  Pressure  Dependence  of  the  Transition  Temperature 

The  development  of  the  energy  gap  in  the  electronic  quasiparticle 

density  of  states  of  the  superconducting  electrode  of  the  tunnel 

junctions  provided  a  very  convenient  method  of  determining  the 

transition  temperature  of  the  alloys  under  investigation.  This 

method  involved  monitoring  the  tunnel  conductance  of  the  samples 

at  zero  bias  (in  the  middle  of  the  energy  gap  region)  and  noting  the 

temperature  at  which  the  sample  conductance  exhibited  a  sharp 

decrease.  Although  the  width  of  the  superconducting  transition 

-2 

has  been  shown  by  Adler  et  al  (1965)  to  be  '5  x  10  K  for  the 
lead-indium  alloy  system  (up  to  26  at  %  In) ,  the  transition 
temperatures  measured  in  this  way  are  still  very  accurate  (better 
than  ±  10  mK) .  This  is  so  since  the  conductance  undergoes  a  very 
large  change  over  this  temperature  and  is  readily  observable  with 
the  high  gain  detection  system  that  was  employed  for  the  whole  of 
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the  tunnelling  experiments  performed. 

The  measurement  of  the  transition  temperature  for  each  sample 

was  repeated  at  least  five  times  and  averaged.  Experimental  accuracy 

of  this  data  was  limited  mostly  by  the  resolution  of  the  electronic 

circuitry  associated  with  the  germanium  resistance  thermometer 

and  to  a  lesser  extent  by  the  uncertainty  in  the  pressure  determination. 

Included  in  Table  1,  which  contains  the  numerical  averages 

obtained,  are  average  values  for  the  ratio  of  room  temperature 

resistance  to  liquid  helium  temperature  resistance  of  the  alloy 

films  (ROQ  /R. )  for  the  respective  compositions.  By  referring 
zy  j  4 

to  figure  11  a  more  precise  value  for  the  average  sample  concentration 
could  be  obtained  than  the  nominal  12  at  %  and  36  at  %  figures  which 
have  been  rounded  to  the  nearest  one  percent  and  which  will  be  used 
henceforth. 

TABLE  1 

Pressure  Dependence  of  the  Transition  Temperature 


Concentration 
(at  %  In  in  Pb) 

R293 

R4 

T  (P=800 
c 

(K) 

p.  s . i. ) 

dT 

c 

dp 

(10"5K  bar  1) 

dlnT 

c 

dp 

(10“6  bar-1) 

12 

3.69 

7.110 

±.008 

-3.19 

±.16 

-4.48 

±.22 

36 

2.18 

6.75 

±.008 

-3.06 

±.15 

-4.53 

±.23 

For  comparison  purposes,  Table  2  lists  the  values  of  T^  determined 
by  other  workers  at  zero  pressure  (for  impurity  concentrations  as 
close  as  possible  to  those  of  interest  here). 
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TABLE  2 

Comparative  Transition  Temperatures  (at  Zero  Pressure) 


Authors 

Concentration 
(at  %  In  in  Pb) 

Tc(P=0) 

(K) 

Adler  et  al  (1967) 

12 

7.097 

Farrel  et  al  (1969) 

10 

7.05 

Farrel  et  al  (1969) 

40 

6.57 

Adler  et  al  (1967) 

40 

6.63 

Wu  (1967) 

40 

6.65* 

^Theoretical 


Several  factors  can  affect  the  critical  temperature  of  a  super¬ 
conductor  when  determined  from  a  thin  film  sample.  The  most  important 
of  these  being  effects  of  structure,  impurity  content  and  stress. 

For  samples  which  are  formed  at  room  temperature  the  crystalline 
structure  is  generally  not  that  different  from  that  of  the  bulk, 
providing  the  film  is  sufficiently  thick,  and  critical  temperatures 
measured  in  this  case  are  indicative  of  those  for  the  bulk  material. 
The  transition  temperature  of  superconductors  deposited  at  low 
temperatures  however  can  exhibit  considerable  enhancement  as  for  the 
case  of  aluminium  for  which  transition  temperatures  up  to  2.5K  have 
been  observed.  The  underlying  reason  is  that  at  low  deposition 
temperatures  the  atom  mobility  of  the  condensate  is  low  and  the 
crystalline  growth  is  inhibited  yielding  an  amorphous  structure. 

An  extreme  example  of  this  effect,  observed  by  Lazaren  et  al  (1957) 
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is  that  of  beryllium  which  in  its  normal  crystalline  state  is  a  non 
superconductor  whilst  if  deposited  at  low  temperature  exhibits  a 
transition  temperature  of  8K.  If  annealed,  the  transition  temperature 
of  a  sample  exhibiting  this  kind  of  Tc  enhancement  reverts  to  that 
of  the  "bulk"  state.  Since  the  samples  used  in  these  experiments 
were  produced  at  room  temperature  and  a  finite  time  necessarily 
elapsed  before  cooling,  structure  enhanced  T^'s  are  not  expected. 

Impurities  constitute  a  more  significant  effect  on  the 
transition  temperature  by  generally  decreasing  the  value  of 
linearly  with  decreasing  electronic  mean  free  path.  According  to 
Mathiessen's  rule  the  decrease  in  in  this  regard  is  linear  and 
given  by:- 

AT  =  Kp 
c  res 

where  K  is  a  constant  associated  with  the  so lute- solvent  pair  and 

p  is  the  residual  resistance.  Since  the  samples  used  were 
res 

produced  by  vacuum  deposition,  the  main  source  of  unwanted  impurities 
would  be  gaseous  in  nature.  Doping  of  indium  with  oxygen  for  example 
depresses  its  transition  temperature.  Depression  of  the  transition 
temperature  due  to  the  inclusion  of  impurities  was  minimised  however 
by  depositing  the  films  from  extremely  pure  source  materials  after 
first  attaining  the  highest  possible  vacuum.  A  liquid  nitrogen  trap 
was  employed  to  assist  in  the  pumping  as  well  as  to  reduce  the 
backstreaming  of  oil  from  the  diffusion  pump. 

Stress  effects  as  a  result  of  differential  contraction  between 
the  substrate  and  the  film  have  the  effect  of  enhancing  the  transition 
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temperature.  This  problem  has  been  investigated  by  Toxen  (1961)  for 
the  case  of  indium  on  fused  quartz  substrates  who  calculated  an 
expression  for  the  change  in  transition  temperature  given  by:- 

AT  =  52  -  750 

C  T  6T 

o 

where  d  is  in  Angstroms.  For  a  film  of  thickness  5000A  this  amounts 

to  an  increase  in  T  of  ~10mK. 

c 

It  is  difficult  to  make  a  quantitative  assessment  of  the  net 

o 

effect  of  the  previous  factors  although  for  thick  films  (>3000A) 

bulk  values  of  T  for  lead  at  least  seem  to  be  in  close  agreement 

c 

with  those  obtained  using  tunnelling  techniques  (see  Franck  &  Keeler 
(1967)). 

It  can  be  seen  from  Tables  1  and  2  that  there  is  good  agreement 
between  the  values  obtained  here  for  the  alloys  with  those  of  other 
workers.  The  last  column  in  Table  1  has  been  included  for  referral 
purposes  in  a  later  section  since  it  provides  a  convenient  form  for 
theoretical  comparison. 

3 .  The  Energy  Gap  and  its  Dependence  on  Pressure 

Determination  of  the  superconducting  energy  gap  was  less  direct 
than  was  the  determination  of  the  transition  temperature.  At  the 
finite  temperatures  employed  here,  the  conductivity  does  not  exhibit 
a  sharp  discontinuity  at  biases  close  to  the  energy  gap  region 
but  appears  as  a  smooth  continuous  function.  Consequently 
a  criterion  has  to  be  established  for  deciding  at  which  point  to 
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evaluate  the  gap.  Figure  12  shows  typical  characteristics  of  a 
Pb^In^  junction  both  at  low  and  high  pressure  recorded  at  a 
temperature  of  ~2.5K  in  which  it  can  be  seen  that  although  a  large 
decrease  in  conductance  occurs,  no  one  obvious  point  is  in  evidence 
for  the  gap  determination. 

One  such  criterion  which  is  widely  used  is  that  proposed  by 
McMillan  &  Rowell  (1969)  for  example  in  which  one  makes  the 
measurement  from  zero  applied  bias  to  the  mid-point  of  the  current 
jump  in  the  current  voltage  characteristic  of  the  superconductor. 

On  the  conductance  curves  that  were  recorded  here  this  point 
corresponds  to  the  inflection  point  on  the  low  energy  side  of 
the  conductance  maximum. 

The  normalised  conductance  of  a  B.C.S.  superconductor  at  finite 
temperatures  has  been  tabulated  by  Bermon  (1964)  as  a  function  of  the 
reduced  parameters  (bias/energy  gap)  and  (energy  gap/kgT)  where  k^ 
is  Boltzman’s  constant  and  T  the  temperature  in  Kelvins.  By  comparing 
the  conductance  plots  obtained  in  these  experiments  with  this  tabulation 
and  identifying  the  one  which  most  closely  approximates  the  experimental 
curve  the  value  of  the  latter  reduced  parameter  was  identified  and 
since  the  temperature  of  the  experimental  determination  was  known, 
it  was  possible  to  calculate  the  gap.  In  employing  this  procedure, 
emphasis  was  placed  on  the  curve  matching  procedure  in  the  higher 
energy  region  where  the  gap  edge  might  be  expected  and  where  the 
error  in  the  conductance  measurements  was  least.  The  criterion  of 
McMillan  and  Rowell  is  implicit  in  this  approach  if,  when  the  parameter 
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FIGURE  12 
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bias/A  assumes  the  value  unity,  the  experimental  and  tabulated 
conductance  values  correspond. 

For  the  example  of  figure  12  this  procedure  yielded  a  value 
of  A  =  1.230meV.  Since  it  would  be  difficult  to  compare  gaps 
measured  at  different  temperatures  with  those  of  other  workers  it 
is  usual  to  convert  these  values  to  their  zero  temperature 
analogues.  This  can  be  done  by  referring  to  the  tabulation  of 
Muhlschlegel  (1959)  which  expresses  the  relationship  between  the 
gap  at  finite  temperature  to  that  at  zero  temperature  again, for  a 
B.C.S.  type  superconductor.  At  the  temperature  of  this  example 
(2.51K)  an  increase  of  1%  in  A  is  predicted  and  we  have,  therefore, 

A  =  1.243meV  for  the  case  at  hand. 

The  whole  of  the  foregoing  has  assumed  the  superconductor  to 
be  of  the  weak  coupling  B.C.S.  type  which  was  not  the  case.  For 
strong  coupling  superconductors  the  energy  gap  is  given  by 
A  =  A^  +  iA2  where  both  A-^  and  are  also  functions  of  energy.  The 
above  method  is,  however,  still  an  excellent  approximation  since 
as  will  be  seen  in  a  later  section,  A2  -  0  and  A^  is  constant  for 
energies  close  to  the  gap  edge. 

The  pressure  induced  shift  in  the  energy  gap  was  obtained  more 
readily  and  with  greater  expected  accuracy  than  the  absolute  value 
of  the  gap  since  for  curves  recorded  at  approximately  the  same 
temperature  the  point  of  comparison  can  be  seen  to  be  non  critical  over 
a  wide  range  of  energy. 

Table  3  summarises  the  results  obtained  for  the  pressure 
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dependence  of  the  energy  gap. 


TABLE  3 

Pressure  Dependence  of  the  Energy  Gap 


Concentration 
(at  %  In  in  Pb) 

A(P=0) 

(meV) 

dA 

dp 

(lCP^meV  bar--*-) 

dlnA 

dp 

(10“6  bar--*-) 

12 

1.32 

-8.14 

o 

• 

+i 

-6.16  ±.30 

36 

1.242 

-8.57 

±.43 

-6.90  ±.36 

Taking  the  logarithm  of  the  B.C.S.  gap  ratio  expression  we  have 


that :- 

din 


2A 


k  T 
B  cJ 


=  dlnA  -  dlnT 
dp 


dp 


dp 


substituting  the  average  values  of  dlnA  and  dlnT  from  Tables  3 

dp  -r— 9- 

dp 

and  1  respectively  in  the  above  expression  gives :- 


din 


/  \ 
2A 


vk_T  , 
v  B  c' 


dp 


and  din 


2A 


KT 

B  c 


=  -  1.68  ±  .52  for  Pbggln^ ^ 


=  -  2.37  ±  .59  for  Pb64In36 


dp 

Both  of  these  results  indicate  a  trend  towards  weaker  coupling 
as  a  function  of  applied  pressure  but  apparently  to  a  lesser  extent 
than  that  found  for  pure  lead  by  Franck  &  Keeler  (1968).  The 


comparitive  rate  of  change  of  the  gap  and  transition  temperature 
can  also  be  estimated  from  the  data  of  Tables  1  and  3  which  give:- 


88 


=  1.375  ±  .53  for  Pbggln^ 

=  1.523  ±  .58  for  Pb,,In0, 

64  36 

Table  4  lists  the  values  of  other  workers  which  are  relevant 
to  those  parameters  determined  here.  In  contrast  to  Table  2y however, 
we  are  now  able  to  compare  the  pressure  dependence  of  our  data  with 
the  very  recent  work  of  Hansen  et  al  (1973)  and  Galkin  et  al  (1973) 


TABLE  4 

Comparative  Pressure  Dependence  of  the  Energy  Gap 


Authors 

Concentration 
(at  %  In  in  Pb) 

A(P=0) 

(meV) 

dA 

dp 

(10”^meV  bar--*-) 

dlnA 

dp 

(10"6  bar"1) 

Adler  et  al 
(1967) 

12 

1.33 

Hansen  et  al 
(1973) 

12 

1.351 

-10.3  ±.06 

-7.63  ±.48 

Galkin  et  al 
(1973) 

12 

1.34 

-8.8  ±.60 

Adler  et  al 
(1967) 

40 

1.21 

Galkin  et  al 

40 

1.22 

-7.9  ±.60 

Hansen  et  al  do  not  appear  to  have  measured  T^  although  in  their 
publication  they  allude  to  a  dlnT^/dp  for  Pbggln.^  of  -5.4  ±2.6  x  10  ^ 
bar  ^  as  being  an  experimental  value.  Using  this  value  (which  pertains  to  Pb) 
with  their  experimentally  determined  dlnA/dp  however  yields  a  value  of 


dlnA 

dlnT 

c 

and  dlnA 
dlnT 

c 
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dlnA/dlnT^  =  1.413  ±  3.1  which,  regardless  of  the  large  tolerance, 

is  in  excellent  agreement  with  the  value  obtained  by  this  author. 

The  decrease  in  the  B.C.S.  gap  ratio  in  strong  coupling 

superconductors  subjected  to  hydrostatic  pressure  (indicative  of  a 

trend  towards  weaker  coupling)  has  been  explained  in  a  series  of 

publications  revolving  around  the  numerical  solution  of  the  Eliashberg 

equations.  The  Eliashberg  equations  (as  was  pointed  out  in  Chapter 

II)  allow  provision  for  the  explicit  form  of  the  phonon  spectrum  of 

a  given  metal  wherein  lies  the  crux  of  the  problem. 

Using  a  Born-von  Karman  force  constant  fit  to  the  dispersion 

curves  obtained  from  inelastic  neutron  scattering  experiments  to 

obtain  information  on  the  lattice  vibrations  and  combining  this 

with  the  Coulomb  pseudopotential  term,  Carbotte  and  Dynes  (1968) 

calculated  the  form  of  the  electron-phonon  interaction  (essential 

to  the  superconducting  state)  and  from  this  calculated  a  form  for 
o 

a  F(oo)  for  use  in  the  Eliashberg  equations  which  led  to  a  value  of 
A  =  1.49  for  lead. 

These  calculations  were  then  extended  to  finite  pressure  by 

2 

Trofimenkoff  and  Carbotte  (1969)  who,  for  the  a  F(oj)  at  finite 

2 

pressure,  simply  scaled  the  zero  pressure  a  F(co)  upward  in  energy 
by  the  observed  average  change  in  sound  velocities.  By  incorporating 

the  experimental  evidence  of  Franck  &  Keeler  for  lead  that  T^  decreases 

h 

according  to  ~ (A) 2  they  were  able  to  predict  qualitatively  the 
observed  reduction  of  Tc  with  pressure. 

Still  within  the  framework  of  the  calculations  of  Carbotte  and 
Dynes,  Trofimenkoff  and  Carbotte  (1970)  again  solved  the  Eliashberg 
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equations  this  time  for  A  and  at  zero  and  at  finite  pressure 

but  with  a  new  scaling  procedure  for  producing  the  pressure 
2 

dependent  a  F (oo)  .  This  scaling,  again  concerned  with  the  electron- 

phonon  interaction,  consisted  of  shifting  of  the  phonon  frequencies 

to  higher  energy  by  an  amount  implied  by  an  average  Gruneisen 

parameter  and  the  calculation  of  a  re-screened  pseudopotential. 

Beside  the  prediction  of  pressure  dependences  of  the  energy  gap  and 

transition  temperature  which  were  in  excellent  agreement  with  the 

experimental  results  Franck  &  Keeler 5 (through  the  relation  dlnA/ 

dlnT^)  they  also  note  that  the  repulsive  Coulomb  part  of  the  electron- 

phonon  interaction  is  fairly  insensitive  to  pressure  and  can  be 

accounted  for  by  a  simple  scaling  factor.  Trofimenkoff  &  Carbotte 

further  predict  a  certain  degree  of  pressure  dependence  in 

dlnA/dlnT  as  can  be  seen  in  Table  5  which  lists  their  values 
c 

together  with  those  obtained  from  experiment. 

TABLE  5 

Relative  Rate  of  Change  of  the  Energy  Gap  and  Transition  Temperature 


Authors 

Concentration 
(at  %  In  in  Pb) 

dlnA 

dlnT 

c 

Trofimenkoff  &  Carbotte 
(1970) 

0 

1.92 

2.18  (P  ~  24.4  kbar) 

Franck  &  Keeler 

0 

2.06  ±  .3 

(1968) 

Hansen  et  al 

12 

1.413  ±  3.1 

(1973) 

This  work 

12 

1.375  ±  .53 

This  work 

36 

1.523  ±  .58 
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There  would  appear  to  be  very  good  agreement  therefore  of 
these  experiments  with  theoretical  predictions.  Carbotte  and 
Vashishta  (1971)  have  extended  these  calculations  to  finite 
temperatures  and  carried  them  out  with  greater  numerical  precision 
although  their  findings  for  the  case  of  strong  coupling  lead 
do  not  differ  significantly  from  those  of  Trof f imenkof f  and  Carbotte. 

The  choice  of  scaling  the  finite  pressure  phonon  spectrum 
to  higher  energy  stems  from  the  work  of  McMillan  (1968)  who  defined 
a  parameter  A  as  a  measure  of  the  coupling  strength  and  showed  that 

A  ^  _1  where  U)  represents  the  frequency  of  the  phonons  responsible 

0)2 

for  the  electron-phonon  coupling. 

The  results  obtained  for  the  alloys  also  indicate  a  decrease 

in  gap  ratio  (and  hence  coupling  strength)  as  a  function  of 

concentration  with  typical  values  of  2A/kTc  being  4.32  and  4.28  for 

the  PbOQIn  and  Pb,.In  alloys  respectively.  This  is  not  too 
oo  1Z  o4  jo 

surprising  since  the  electron/atom  ratio  has  decreased  in  each  case 
resulting  in  a  reduction  in  the  screening  of  the  ion  sites  and  thus 
a  corresponding  reduction  in  the  strength  of  the  electron-phonon 
interaction. 

Further  effects  on  the  gap,  T^  and  the  gap  ratio  include  those 
effects  mentioned  in  the  previous  section  (for  T  )  and  possible  gap 
anisotropy  effects.  The  transition  temperature  can  also  be  dependent 
on  crystal  alignment  in  the  film.  The  net  effect  of  these  latter 
factors  is,  for  all  intents  and  purposes,  inestimable. 


. 


'  ■ 
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A .  Pressure  Dependence  of  the  Phonon  Spectrum 

A  great  deal  of  information  regarding  the  phonon  spectra  of 
the  alloys  can  be  obtained  from  the  tunnel  conductance  characteristics 
some  of  which,  namely  the  location  and  pressure  shifts  of  the 
principle  phonon  modes  (and  hence  the  associated  Griineisen  parameters), 
is  observable  in  a  very  direct  manner. 

Further  information  can  be  obtained  somewhat  less  directly  in 
a  procedure  (to  be  described)  which  involves  the  numerical  solution 
of  the  Eliashberg  equations.  This  latter  approach  involves 
extensive  computation  but  produces  as  a  result  an  overall  view  of  the 
phonon  spectrum,  the  renormalization  function  and  the  gap  function 
as  well  as  such  important  parameters  as  A,  (the  coupling  strength) 
and  p*  (the  Coulomb  pseudopotential  term) . 

The  two  different  approaches  employed  in  dealing  with  the 
data  of  this  section,  therefore,  suggest  a  natural  subdivision  for 
its  analysis. 

a.  Directly  Observed  Principle  Phonon  Characteristics. 

It  was  shown  in  Chapter  II  that  the  normalised  tunnel  conductance 
was  in  one  to  one  correspondence  with  the  electronic  quasiparticle 
density  of  states  and  further  that  structure  in  the  density  of  states 
was  a  reflection  of  the  phonon  characteristics  of  the  metal  in 
question.  Since  phonon  structure  is  measured  from  the  gap  edge, 
the  tunnel  characteristics  and  their  derivatives  appearing  henceforth 
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are  displayed  so  that  their  gap  edge  energies  coincide.  In  this  way, 
shifts  in  phonon  energies  are  directly  observable.  Figures  13  and  14 
show  characteristics  for  typical  Pbnr>Inno  and  Pb^.In-.  tunnel  junctions 
together  with  their  B.C.S.  weak  coupling  counterparts  (which  do  not 
contain  the  specific  nature  of  the  phonon  spectrum).  The  general 
decrease  in  the  sharpness  of  the  structure  of  these  curves  with 
increasing  indium  concentration  is  consistent  with  the  trend  towards 
weaker  coupling  that  was  observed  previously  in  the  gap  ratio  at 
zero  pressure.  Critical  points  in  the  phonon  spectrum  corresponding 
to  the  predominant  phonon  modes,  manifest  themselves  in  the  conductance 
curves  as  points  of  inflection  which  can  be  easily  resolved  by 
differentiation  whereupon  each  point  of  inflection  gives  rise  to 
a  peak. 

The  phonon  structure  observed  in  this  way  was  found  to  be 
sharper  with  decreasing  temperature  and  thus  the  lowest  temperatures 
that  could  be  maintained  fairly  constant  were  employed  in  these 
determinations.  With  decreasing  temperature  however  the  possibility 
of  instrumental  smearing  of  the  data  becomes  more  prevalent.  A  temperature 

of  1.4K  for  example ,  corresponds  to  a  thermal  smearing  of  ~120yV  of 
the  Fermi  surface  and  thus  to  ensure  that  the  smearing  of  data  was 
thermally  rather  than  instrumentally  limited,  applied  modulation 
signals  should  be  as  small  as  possible  and  certainly  less  than  that 
corresponding  to  the  temperature  at  which  the  measurements  are  to  be 
made.  For  the  second  derivative  measurements  here  (and  also  for  the 
gap  determinations)  applied  modulation  signals  were  of  the  order  of 
lOyV  -  20yV  r.m.s.  It  is  note  worthy  that  during  the  measurement 
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Comparison  with  the  B.C.S.  Theory  (P=0) 
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FIGURE  14 
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of  the  energy  gaps,  it  was  observed  by  this  author  that  the 
conductance  maximum  (at  least)  was  definitely  dependent  on  the 
magnitude  of  the  applied  modulation  signal  for  signals  in  the 
range  >100yV  r.m.s.  At  the  particular  point  of  observation  the 
shift  was  observed  to  be  of  the  order  .4%  for  an  increase  in 
modulation  of  70yV  r.m.s.  above  the  lOOyV  level. 

Figure  15  shows  a  typical  conductance  characteristic  together 

with  the  associated  second  derivative  in  this  case  for  a  Pb..,In~. 

64  36 

sample  as  an  illustration  of  the  general  method  of  locating  the 
principle  phonon  peaks. 

Figures  16  and  17  show  the  normalised  conductance  characteristics 
for  typical  samples  of  Pbggln.^  and  Pb^In^  that  were  obtained  as 
a  function  of  pressure.  By  comparing  qualitatively  the  figure  17 
with  the  figure  16  we  again  see  the  general  washing  out  of  the 
phonon  structure  with  increasing  indium  content.  Less  obvious, 
but  still  apparent,  by  examining  figures  16  and  17,  is  the  general 
shift  in  the  phonon  spectrum  to  higher  energies.  The  structure 
amplitude  reduction  which  occurs  as  a  function  of  pressure  (as 
opposed  to  indium  concentration)  is  not  easily  observable  from 
these  curves  as  a  consequence  of  the  scale  used. 

The  principle  difference  in  the  phonon  spectrum  of  the  lead- 
indium  alloy  system  from  that  of  pure  lead  is  the  existence  of  the 
localised  mode  of  vibration.  The  power  of  the  method  of  differentiating 
the  conductance  curve  to  locate  the  phonon  structure  is  emphasised  in 
figure  18  which  shows  the  high  energy  region  for  a  typical  Pb^^In^g 
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First  and  Second  Derivatives  for  a  Pb  In  Sample  (P=0) 

“  1  . . .  ■  i  ,  64  36  ~ 
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FIGURE  18 
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sample.  In  this  figure  the  shifts  in  the  longitudinal  and  localised 
modes  are  clearly  resolved  as  a  function  of  pressure. 

In  resolving  the  predominant  phonon  peaks  a  difficulty  is 
encountered  in  the  low  energy  region  where  two  transverse  phonon 
branches  exist  at  slightly  different  energies.  This  has  the  overall 
effect  of  broadening  the  peak  while  also  adding  slight  structure. 

By  using  the  minimum  in  the  peak  obtained  for  the  transverse  mode 
(as  was  done  for  the  longitudinal  and  localised  modes)  what  has  in 
fact  been  measured  is  the  higher  energy  branch. 

Table  6  summarises  the  results  pertaining  to  the  phonon  energies 
and  their  pressure  dependence  in  which  the  subscripts  t,  1  and  i  refer 
to  the  transverse,  longitudinal  and  localised  (impurity)  modes 
respectively. 


TABLE  6 


Pressure  Dependence  of  the  Principle  Phonon  Peaks 


Concentration 

Peak 

Location  (P=0) 

Pressure  Dependence 
(10“^meV  bar“l) 

(at  %  In  in  Pb) 

(meV) 

dco. 

dU)^ 

dp 

dto . 

l 

dp 

“i 

dp 

12 

4.30 

8.65 

9.86 

2.74  ±.69 

5.21  ±.52 

4.33  ±.65 

36 

3.67 

8.79 

10.43 

2.24  ±.56 

4.97  ±.50 

4.01  ±.6 

Hansen  et  al  (1973)  and  Galkin  et  al  (1973)  have  very  recently  obtained 
results  on  the  phonon  shifts  as  a  function  of  pressure  for  the  lead- 
indium  alloy  systems  and  their  results  are  summarised  in  Table  7. 


Comparative  Pressure  Dependence  of  the  Principle  Phonon  Peaks 
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By  reaching  lower  temperatures  (at  which  the  aluminium  electrode 
was  superconducting)  Galkin  et  al  have  been  able  to  resolve  the  two 
transverse  modes  through  increased  sensitivity.  Their  results  were 
recorded  for  pressures  between  9kbar  and  llkbar  which  suggests  (using 
the  data  cited  by  Trofimenkoff  and  Carbotte  (1969))  that  the  tempera tur 
of  their  investigations  was  in  the  neighbourhood  of  ~0.7K.  The 
general  agreement  between  the  results  obtained  here  and  those  of 
Hansen  et  al  and  Galkin  et  al  is  excellent. 

Attention  is  now  turned  towards  the  calculation  from  Table  6 
•  • 

of  the  relevant  Gruneisen  constants.  The  Gruneisen  rule  states 
that  at  low  temperatures :- 


Yg  =  Va 


KC 


V 


6-1 


where  V  is  the  volume,  a  the  thermal  expansion  coefficient,  the 
specific  heat  at  constant  volume  and  K  the  compressibility.  If  we 
consider  this  in  terms  of  the  normal  modes  of  vibration  of  the 
substance  it  can  be  expressed  as:- 


Y  = 


Zy.C  . 
.  1 1  vx 
i 

SC  . 

.  vi 
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where  y^  are  the  individual  mode  gammas  defined  by:- 


y.  =  -  dlnoo .  = 
3  _ i 

dlnV 


ffiohl  B 

l 


l6P 


6-3 


03. 


where  the  ol  are  now  the  frequencies  associated  with  the  normal 

modes  of  vibration  and  B  is  the  bulk  modulus  defined  by  B  =  -  dp/dlnV, 
•  • 

The  Gruneisen  constants  for  the  predominant  phonon  modes  can 
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therefore  be  determined  and  are  given  in  Table  8. 

TABLE  8 

Gruneisen  Gammas  of  the  Predominant  Phonon  Modes 


Concentration 

Gruneisen  mode 

y's 

(at  %  In  in  Pb) 

rt 

A 

A 

12 

3.11  ±.78 

2.94  ±.29 

2.14  ±.32 

36 

2.98  ±.75 

2.76  ±.28 

1.88  ±.28 

The  data  of  Bridgman  (1954)  obtained  at  room  temperature  shows 
that  the  bulk  modulus  of  Pb^I^^  is  not  significantly  different  from 
that  of  pure  lead  and  thus,  in  the  absence  of  any  data  for  the  alloy 
at  low  temperatures,  the  Gruneisen  gammas  have  been  determined  here 
using  a  value  of  4.88  x  10~*  bar  for  the  bulk  modulus  of  pure  lead  at 
4K  as  determined  by  Waldorf  and  Alers  (1961). 

Listed  in  Table  9  are  the  values  for  the  Gruneisen  y's  as 
calculated  from  the  data  of  Galkin  et  al. 


TABLE  9 

Gruneisen  Gammas  Calculated  from  the  Data  of  Galkin  et  al  (1973) 


Concentration 

Gruneisen  mode 

Y*  s 

(at  %  In  in  Pb) 

Yt 

^1 

Yi 

12 

3.90  ±.49 

3.03  ±.24 

2.73  ±.25 

40 

3.66  ±.48 

2.73  ±.34 

1.32  ±.34 
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Agreement  of  the  results  obtained  here  is  also  very  good  and 
always  within  the  experimental  error.  It  is  noteworthy  that  for  the 
40  at  %  sample  of  Galkin  he  no  longer  was  able  to  resolve  two 
transverse  modes  and  that  the  structure  he  observed  in  the  second 
derivative  became  unduly  washed  out  which  leads  this  author  to 
believe  that  perhaps  the  concentration  was  higher  than  that 
assigned . 


b .  Inversion  of  the  Eliashberg  Gap  Equations. 

2 

By  guessing  the  parameters  a  (co )  and  F(oo),  relating  to  the 
phonon  spectrum  of  the  metal,  it  is  in  principle  possible  to  solve 
the  Eliashberg  coupled  integral  equations  described  in  Chapter  II 
for  and  Z(oo)  and  hence  obtain  the  gap  function  A  (co)  through 

the  relationship :- 

A(co)  =  ^  (oo)  / Z  (co)  2-66 

Subsequently  one  can  use  this  to  calculate  the  normalised  density 
of  states  for  a  superconductor  using  the  relationship : - 

D  (E)  =  Re  |  eV  |  2-74 

D(ef)  [(eV)2  -  k2(ui)]% 

Since  the  normalised  density  of  states  was  earlier  shown  to  be  in 
1:1  correspondence  with  the  normalised  tunnel  conductance,  the  results 
of  such  a  calculation  can  be  compared  directly  with  experiment. 

This  procedure  was  carried  out  by  Rowell  and  Kopf  (1965)  and  others 
with  excellent  results. 

McMillan  (1965  and  1969)  and  Rowell  however  have  devised  an 
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innovative  approach  to  solving  the  Eliashberg  equations  (and  also 
using  equation  2-74)  which  uses  experimentally  obtained  input  data 
for  the  energy  gap  and  density  of  states  and  yields  (after  considerable 
computation)  the  form  of  the  phonon  spectrum  and  associated  parameters. 
Because  the  aim  of  this  procedure  is  to  gain  information  regarding 
the  phonon  spectrum  rather  than  the  density  of  states  (as  was  the 
previous  method  of  solution)  it  has  come  to  be  known  as  an  inversion 
process . 

Briefly  the  method  is  as  follows.  A  zero  order  guess  is  first 

2 

made  for  the  function  a  (uj)F(oij),  (the  Coulomb  pseudopotential)  and 

A.  The  first  two  of  these  parameters  are  based  on  intelligent  insight 

while  the  initial  input  for  A  is  the  gap  at  the  gap  edge  as  obtained 

from  experiment.  On  the  basis  of  this  information,  the  gap  equations 

are  first  solved  by  iteration  for  A(cu)  from  which  the  density  of  states 

is  calculated.  This  calculated  density  of  states  is  then  compared 

with  that  obtained  by  experiment  and,  depending  on  the  difference 

2 

between  the  two,  a  revised  guess  is  made  for  a  (oo)F(oj).  This 

procedure  is  repeated  as  many  times  as  is  necessary  until  the  agreement 

between  the  calculated  and  experimental  densities  of  states  is 

2 

acceptable  at  which  time  the  best  form  of  a  (co)F(to)  will  have  been 
arrived  at.  Between  successive  iterations  for  A(co),  the  parameter 
U  is  adjusted  so  that  the  calculated  gap  function  (at  the  gap  edge) 
is  consistent  with  the  experimental  value. 

Output  from  this  kind  of  analysis  includes  the  form  of  the 
Coulomb  pseudopotential  y*  (note  y*  =  N(0)Uc)  as  well  as  of  course 
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a  (oj ) F (jLO )  (from  which  the  electron  phonon  coupling  strength  \  and 
various  average  phonon  energies  can  be  derived). 

Such  an  analysis  has  been  carried  out  for  the  Pb^In^  system. 
As  input  data  to  the  computer  programme  the  function  (a/a^  _  c  )-l 
was  required  where  the  a  refer  to  normalised  conductances  and  in 
particular 


a 


V 


B .  C .  S  . 


(V2  - 


>2)* 


In  obtaining  a  for  the  alloy  from  (a^g/cs^),  it  was  first  necessary 
to  apply  a  correction  to  the  absolute  value  of  before  obtaining 

the  low  pressure  data.  This  arose  from  the  fact  that  and  o.n, 

were  determined  at  slightly  different  pressures.  At  the  present 
time  it  is  sufficient  to  say  that  this  correction  is  entirely 
justified  and  will  be  dealt  with  in  greater  detail  in  the  following 
section  (5). 

Figure  19  shows  the  input  data  used  for  the  inversion  programme 
while  figures  20,  21  and  22  show  the  resulting  phonon  spectrum,  gap 
function  and  re-normalisation  function  respectively. 

To  the  knowledge  of  this  author,  no  other  pressure  dependent  data 
is  available  with  which  to  compare  the  results  at  this  concentration. 
Wu  (1967)  has  obtained  the  phonon  spectrum  at  zero  pressure  for  a 
Pb^In^Q  alloy  which  compares  favourably  with  that  of  figure  20 
in  a  qualitative  sense  although  the  localised  mode  appears  to  be 
sharper  than  might  be  expected  (Wu's  data).  The  gap  function  and 
re-normalisation  function  that  have  been  calculated  here  show 
considerably  less  structure  than  those  for  lower  concentration 
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Energy  (meV) 


a2(co)F(co) 
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FIGURE  20 


FIGURE  21 


110 


Energy  ( meV ) 


FIGURE  22 
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together  with  a  general  slight  decrease  in  magnitude  indicative  of 
the  trend  to  weaker  coupling  as  a  function  of  concentration.  Also 
in  evidence  is  the  trend  to  weaker  coupling  as  a  function  of  pressure. 

As  can  be  seen  from  figure  20,  the  overall  phonon  structure  has 
decreased  in  magnitude  and  definitiveness.  This  is  especially 
noticeable  in  the  case  of  the  localised  mode.  A  drastic  decrease 
in  the  amplitude  of  this  peak  is  evident  however  as  the  indium 
concentration  is  increased  from  ~12  at  %  to  ~40  at  %  as  can  also 
be  seen  in  the  data  of  Galkin  et  al .  Sood  (1972)  has  investigated 
the  effect  of  alloying  on  the  phonon  spectrum  of  anc^ 

concludes  that  the  localised  impurity  band  associated  with  the 
impurity  atoms  of  indium  consists  in  fact  of  two  peaks.  The  first 
of  these  peaks  (lower  energy)  he  finds  to  occur  at  constant  energy 
(9.57  ±.03  meV)  and  with  width  that  varies  linearly  with  concentration 
while  the  second  peak  not  only  varies  in  width  with  concentration  but 
also  shifts  to  higher  energies  as  a  function  of  increasing  indium 
concentration  (for  concentration  >2  at  %  In).  The  first  of  these 
two  peaks  he  associates  with  the  vibration  of  single  substitutional 
impurity  indium  atoms  in  the  lead,  lattice  while  the  second  peak  is 
believed  to  be  due  to  the  vibration  of  pairs  of  indium  atoms. 

This  effect  certainly  would  be  inherent  in  the  lack  of 
definition  of  the  localised  mode  but  to  an  even  greater  extent  since 
not  only  is  it  a  poor  approximation  to  say  that  the  indium  atoms 
are  localised  but  it  becomes  difficult  to  evaluate  the  effect 
of  the  change  in  the  force  constants  at  such  a  concentration.  In 
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short  it  is  somewhat  surprising  that  the  localised  mode  is  observed 
at  all. 

McMillan  (1968)  defines  a  dimensionless  electron  phonon  coupling 
constant  as:- 


A  =  2 


a  (to) F  (to)  do) 


6-4 


03 


and  within  the  framework  of  the  Eliashberg  strong  coupling  gap 

equations  derives  a  formula  for  the  transition  temperature  for  a 

strong  coupling  superconductor  given  by:- 

T  =  <to>  exp  [-1.04(1  +  A)/  (A  -  y*  -.62Ay*)]  6-5 

C  1.20 


From  equation  6-5  the  following  expression  for  the  pressure  dependence 
of  the  transition  temperature  can  be  derived :- 


dlnT  =  dln<03>  +1.04 

c 


dp 


dp 


1  +  .38y* 

(A-y*-.62Ay*)2J 


dA 

dp 


where  <U3>  is  an  average  phonon  frequency  defined  by:- 


<03>  = 


a  (oj)F (qj)dqj 

a  (w) F(oj)daj 


03 


For  convenience  we  introduce  the  parameter  A  - 

i 

equation  6-7  can  be  expressed  as:- 
2 

<03>  =  2A 


6-6 


6-7 


a  (03 ) F  (o) )  and  thus 


6-8 


McMillan  suggests  that  for  f.c.c.  lattices  a  good  approximation  to 
<U3>  is  an  average  of  the  transverse  and  longitudinal  phonon  frequencies 

i. e.  :- 
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<C0>  -  0)  =  +  OJ^) 


6-9 


differentiating  equation  6-9  with  respect  to  pressure  gives :- 


dln<oo>  = 

1  )  ^t  + 

r  i  ) 

dp 

k+uiJ  dP 

« — i 
3 

+ 

+j 

3 

doo. 
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Table  10  summarises  the  data  that  was  obtained  using  the  inversion 
technique  for  the  Pb^In^  alloy. 

Using  the  experimental  data  of  Table  6  in  conjunction  with 
equation  6-10  to  obtain  the  first  term  in  equation  6-6  and  the  data 
of  Table  10  for  the  second  term  yields  a  value  of:- 

6-11 


dlnT  =  -3.01  x  10  6  bar  1 

_ c 

dp 


This  compares  with  the  direct  experimentally  determined  value  of:- 

dlnT  =  -4.53  x  10"6  bar"1  6-12 

c 


dp 

If  instead  we  use  the  expression  6-8  we  have  that:- 

2 

dln<qj>  =  dlnA  -  dlnA 
dp  dp  dp 

using  this  expression  (with  values  from  Table  10)  in  the  McMillan 
equation  6-6  gives :- 

dlnT  =  -7.18  x  10~6  bar"1  6-13 

_ c 

dp 

The  agreement  obtained  between  the  experimentally  determined 
dlnT^/dp  and  that  calculated  using  McMillan's  formula  is  much  better 
in  the  first  case  (equation  6-11)  than  in  the  second  case  (equation 
6-13),  the  difference  revolving  around  the  definition  of  <oo>.  In 
this  regard  it  is  interesting  to  note  that  better  agreement  was 
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obtained  for  the  value  which  was  less  dependent  on  the  "inversion 
acquired"  data. 

Although  the  McMillan  equation  for  is  in  itself  approximate, 
this  author  attributes  the  major  proportion  of  the  disagreement  to 
the  inversion  procedure. 

The  experimentally  determined  density  of  states  used  as 

input  to  the  inversion  programme  is  particularly  inaccurate  in  the 

low  energy  region  and  thus  it  is  usual  to  make  some  kind  of  assumption 

here.  The  basis  for  this  assumption  has  long  been  the  subject  of 

2 

some  conjecture  as  to  whether  a  (aj)F(a))  depends  linearly  or 

quadratically  on  oo.  A  very  recent  paper  by  Allen  and  Silberglitt 

2 

(1974)  suggests  that  at  low  frequency,  a  ((jo)F(oo)  is  linear  in  (0. 

2 

Such  factors  will  clearly  affect  the  computed  a  (to)F(oo)  and  y* 
upon  which  the  parameters  of  Table  10  and  calculations  employing 
them  depend . 
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5 .  The  Tunnel  Junction  as  a  Pressure  Sensing  Device 

For  the  entirity  of  this  research  work,  the  pressure  that  was 
ultimately  applied  to  the  samples  was  determined  by  the  sample 
itself  as  alluded  to  earlier.  So  reliable  and  reproducible  were 
the  pressure  responses  of  the  tunnel  junctions  that  their  use  as 
pressure  gauges  in  their  own  right  merits  further  attention. 

In  dealing  with  the  theory  of  tunnelling  (Chapter  II  section  6) 
it  was  shown  that  for  the  normal-normal  metal  case  that  the  net 
tunnel  current  which  flowed  as  a  result  of  an  applied  potential  V 
was  given  by:- 

I  a  D2T2eV  6-14 

2 

where  D  represents  the  product  of  the  electron  density  of  states 
functions  (close  to  the  Fermi  surface) ,  e  the  electronic  charge 
and  T  some  average  matrix  element  describing  the  probability  that  an 
electron  will  in  fact  tunnel  from  one  metal  to  the  other. 

By  considering  the  simple  case  of  one  dimensional  tunnelling,  one 
can  obtain  an  estimate  of  the  quantity  T  by  solving  the  Schrodinger 
equation  for  a  simple  potential  barrier  to  find  that:- 

T  a  exp-[  (2m{U-E})  2t/ti]  6-15 

where  E  is  the  kinetic  energy  of  the  electron,  U  is  the  height  of  the 
barrier  and  U>E;  t  represents  the  thickness  of  the  barrier.  Combining 
equations  6-14  and  6-15  we  get:- 

I  cx  OjVexpf-^  2t]  6-16 

o 

where  ^  the  work  function  in  eV  replaces  (U-E) ,  t  is  in  A  and  C  are 
constants.  Using  Ohm's  law  the  tunnel  resistance  is  then  given  by:- 
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6-18 


6-17 


By  applying  pressure  to  the  tunnel  junctions  we  have  been  able  to 
test  the  validity  of  this  relationship  since  a  reduction  in  t  should 
produce  an  exponential  increase  in  P. 

Results  of  this  investigation  are  shown  in  figures  23  and  24 
for  typical  samples  of  each  composition.  In  each  case  the  data  has 
been  plotted  linearly  and  logarithmically  to  illustrate  the 
exponential  behaviour. 

Since  the  circuitry  employed  in  the  determination  of  the  tunnel 
resistance  was  capable  of  resolving  resistance  changes  of  a  few  parts 
in  10^,  experimental  uncertainty  of  this  data  lies  in  the  pressure 
determinations.  The  smallest  graduation  of  the  Heise  Bourdon  gauge 
was  100  p.s.i.  while  the  manganin  gauge  supplied  with  the  compressor 
was  accurate  to  better  than  1%.  The  uncertainty  of  the  pressure  values 
is  therefore  estimated  to  be  less  than  1%. 

It  can  be  clearly  seen  from  figures  23  and  24  that:- 


6-19 


and  hence  using  equation  6-18  that:- 


-P  oc  C6^  2t  +  C? 


6-20 


Equation  6-19  is  precisely  what  one  would  expect  assuming  t  to 
change  linearly  with  pressure. 

For  the  most  part  the  tunnel  junctions  were  used  to  determine 


the  pressure  in  the  bomb  in  the  ~1K  —  9K  temperature  range  although 
their  calibration  (against  the  two  external  gauges)  took  place  at 
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FIGURE  23 


FIGURE  24 
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~30K.  In  view  of  this,  the  dependence  of  the  zero  bias  tunnel 
resistance  on  temperature  was  checked  by  recording  data  at  ~7K  and 
also  at  ~77K.  Variation  of  R^  over  this  range  was  found  to  be  from 
.2/  to  .5/  and  negligible  in  comparison  to  the  ~30%  change  invoked 
by  a  pressure  change  of  "''50,000  p.s.i. 

Reproducibility  of  the  data  when  subjected  to  pressure  cycling 
was  found  to  be  excellent,  the  tolerance  after  successive  runs  being 
better  than  ±  .5%.  When  maintained  at  temperatures  <77K  the 
characteristics  of  the  tunnel  junctions  tested  remained  unchanged  for 
periods  in  excess  of  two  weeks  and  although  prolonged  exposure  of  the 
tunnel  junction  to  room  temperature  caused  deterioration,  the  only 
outcome  of  a  brief  exposure  was  the  need  for  re-calibration. 

Beside  its  most  widely  acknowledged  use  therefore  as  an  analytical 
tool,  the  tunnel  junction  also  exhibits  great  potential  as  a  sensitive 
pressure  monitoring  device  with  application  over  a  wide  range  of 
temperatures  and  in  both  solid  and  fluid  media.  Throughout  the  course 
of  this  work  the  tunnel  junction  itself  has  proved  an  invaluable  tool 
since  having  filled  the  sample  chamber  with  solid  helium  four  -  and 
thereby  rendering  the  external  pressure  gauges  ineffective  -  it 
provided  the  only  convenient  method  of  monitoring  further  pressure 
changes  (at  constant  volume)  which  might  occur  as  a  result  of 
temperature  changes  for  example.  The  following  example  taken  from 
this  work  illustrates  not  only  the  sensitivity  of  the  tunnel  junction 
but  its  agreement  with  thermodynamical  considerations  of  the 
pressure  transmitting  medium. 
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In  order  to  quench  superconductivity  within  the  lead-indium 
alloy  electrode  of  the  tunnel  junction  and  thereby  obtain  the  normal- 
normal  tunnel  resistance  as  a  function  of  d.c.  bias,  the  temperature 
of  the  sample  chamber  was  raised  to  ~9K  whereas  the  superconducting 
normal  characteristic  was  recorded  at  ~1.5K.  Theoretical  predictions 
indicate  that  at  constant  pressure  and  for  energies  well  beyond  the 
superconducting  gap  region  these  characteristics  should  be  coincident, 
however,  in  view  of  the  temperature  difference  between  these  two 
measurements  (and  accompanying  constant  volume  pressure  change) 
one  might  expect  to  observe  a  certain  degree  of  non-coincidence. 

Figures  25  and  26  show  typical  examples  of  these  characteristics 
for  a  Pb^  In^6  sample  both  at  high  pressure  (~50,000  p.s.i.)  and  at 
low  pressure  (~800  p.s.i.).  Despite  the  temperature  difference  in 
the  measurements,  it  can  be  seen  that  for  the  high  pressure  case 
the  characteristics  are  in  fact  coincident  suggesting  that  a  negligible 
pressure  change  has  taken  place.  The  helium  four  melting  curve 
below  taken  from  Spain  &  Segall  (1970)  shows  that  a  solid  formed 
at  50,000  p.s.i.  (~3.4  K  bar)  remains  in  the  solid  state  at  9K 
and  also  that  in  general,  little  change  in  pressure  occurs  in  the 
solid  as  a  function  of  temperature. 
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FIGURE  25 


Resistance  ( ohm)& 
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FIGURE  26 
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T  1  4 

Isochores  of  Solid  He  at  High  Pressures 
(figures  on  curves  are  molar  volumes  in  cm^) 


Referring  to  the  more  detailed  results  of  table  11  ,  also  taken 
from  Spain  &  Segall,  it  is  indicated  for  example  that  in  a  solid 
formed  at  3441  bar,  a  pressure  change  of  the  order  of  .2%  occurs  as 
a  result  of  warming  the  sample  the  necessary  9K  to  achieve  the  normal- 
normal  state.  Thus  we  observe  exactly  what  would  be  expected. 

For  the  low  pressure  measurements  however,  which  were  carried 
out  at  800  p.s.i.,  the  9K  rise  in  temperature  causes  a  change  of 
state  of  the  pressure  transmitter  from  a  solid  to  a  fluid  when 
recording  the  normal-normal  characteristics  as  can  be  seen  from 
the  graph  below  taken  from  Edwards  and  Pandorf  (1965):- 
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Isochores  of  Solid  He^  at  Low  Pressures 
(figures  on  curves  are  molar  volumes  in  cm^) 


To  estimate  the  pressure  in  the  cell  therefore  at  ~9K  it  is  necessary  to 
refer  to  isochores  in  fluid  helium.  The  data  shown  in  table  12 
corresponds  to  such  isochores  which  have  been  taken  from  Hill  and 
Lousnasmaa  (1959).  A  solid  formed  at  ~800  p.s.i.  corresponds  to  a 
molar  volume  of  19.5  cm^  and  to  obtain  the  pressure  in  the  fluid 
at  this  molar  volume,  a  linear  extrapolation  (extrapolated  values 
appear  in  parentheses)  has  been  performed  on  the  above  data  which 
shows  the  pressure  to  be  ~1900  p.s.i.  i.e.  an  increase  of  ~1100  p.s.i. 

Employing  standard  curve  fitting  techniques  to  the  experimental 
data  of  figure  23 s  the  expected  resistance  difference  between  the 
normal  and  superconducting  characteristics  which  corresponds  to  this 
pressure  change  is  ~.7  ohm  and  as  can  be  seen  in  figure  26  this  is 
exactly  what  was  observed.  Thus  we  can  see  that  there  is  excellent 
agreement  between  the  pressure  determined  by  the  tunnel  junction  and  that 
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TABLE  12 

Pressure  Values  (in  p.s.i.)  Along  Isochores  of  Fluid  Helium 


N\.  Molar  Vol 

N\s^cin3) 

21.55 

20.78 

20.06 

19.5 

19.39 

Temp 

(K)  \ 

3.0 

616 . 1 

801.9 

1021.7 

3.25 

635.8 

821.0 

1041.5 

(1257.9) 

1312.0 

3 . 5 

654.9 

840.8 

1062.8 

(1279.8) 

1334.0 

3.75 

674.7 

861.4 

1084.9 

(1302.4) 

1356.8 

4.0 

695.3 

882.7 

1107.6 

(1326.4) 

1381.1 

4.25 

716.6 

905.5 

1131.3 

(1351.1) 

1406.1 

4.50 

738.7 

929.0 

1155.4 

(1376.5) 

1431.8 

4.75 

761.5 

953.3 

1180.4 

(1402.7) 

1458.2 

5.0 

785.0 

978.3 

1206.1 

(1429.6) 

1485.4 

5.25 

809.2 

1003.3 

1232.6 

(1457.2) 

5.5 

834.2 

1029.0 

1259.8 

(1485.6) 

6.0 

884.9 

1082.7 

1315.7 

(1542.6) 

6.5 

938.6 

1138.5 

1373.7 

(1601.9) 

7.0 

993.7 

1195.8 

1433.3 

(1662.6) 

7.5 

1050.3 

1254.6 

1494.3 

(1724.8) 

8.0 

1107.6 

1314.2 

(1556.0) 

(1787.7) 

9.0 

1225.2 

1436.9 

(1680.9) 

(1913.8) 

. 
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predicted  by  the  thermodynamics  of  the  pressure  transmitter.  Using 
the  junction  in  this  fashion  it  has  been  possible  to  test  whether 
the  constant  pressure  freezing  process  used  did  in  fact  produce 
hydrostatic  pressure  in  the  sample  chamber. 

Although  the  previous  thermodynamical  argument  verifies  in 
general  that  the  pressure  was  indeed  hydrostatic,  a  more  quantitative 
estimate  was  obtained  by  producing  two  junctions  on  the  same  substrate 
a  fixed  distance  apart  and  monitoring  the  pressure  at  these  points. 

By  so  doing  the  pressure  was  found  to  vary  by  <1%  over  a  distance 
of  ~1.5  cm  after  constant  pressure  freezing. 

With  particular  reference  to  this  work,  the  inescapable 
conclusion  is  that  the  tunnel  junction  is  not  only  an  accurate 
method  of  determining  the  pressure  in  the  solid  and  in  the  fluid 
but  it  does  so  identically  at  the  point  of  interest. 
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Pressure  Dependence  of  the  Density  of  States 


For  a  strong  coupling  superconductor,  the  electronic  density 
of  states  at  the  Fermi  surface  is  renormalised  according  to:- 

Nf  =  Nb  (1  +  X)  6-21 


where  X  is  the  coupling  strength  and  1  +  X  =  Z,  the  renormalisation 
constant,  and  Z  =  m*/m  where  m*  and  m  are  the  effective  mass  and 


mass  of  the  electrons  respectively.  N  denotes 
density  of  states. 


tne  Dana  structure 


The  pressure  dependence  of  the  density  of  states  therefore  is, 
using  equation  6-21:- 


dlnN„  =  dlnN  +  dlnZ 
t  a  — — 

dp 


6-22 


dp  dp 

The  pressure  dependence  of  X  for  36  at  %  In  in  Pb  was  determined 

in  these  experiments  from  which  (using  the  relationship  (l/x)dx  = 

dp 


dlnx/dp  and  the  fact  that  Z  =  1  +  X) 

dlnZ  =  -4.32  x  10  ^  bar  ^ 
dp 

or  since  -1  dV  =  -dlnV  =  K 
V  dp  dp 

dlnZ  =2.11 
dlnV 


6-23 

6-24 

6-25 


where  is  the  isothermal  compressibility  of  bulk  lead.  The 
quantity  dlnN^/dp  has  been  measured  by  several  workers  which  are 
shown  in  Table  13.  Also  appearing  in  this  table  are  the  values  for 
the  pressure  dependence  of  the  band  structure  density  of  states 
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using  equation  6-22  and  the  value  of  dlnZ/dp  obtained  from  these 
experiments.  It  is  expected  that  dlnZ/dp  for  the  alloy  does  not 

differ  considerably  from  that  of  pure  lead  to  which  the  values  of 
column  3  pertain. 


TABLE  13 

Pressure  Dependence  of  the  Density  States  (Per  Unit  Volume) 


Author 

Experimental 

Method 

dlnN„ 

r 

,dp 

(1CT6  bar“l) 

dlnN^ 

D 

dp 

(10-6  bar“d) 

Andres  (1961) 

critical  field 

+.62  ±  3.08 

4.94 

White  (1961) 

thermal  expansion 

-1.44  ±  1.02 

2.88 

Andres  (1964) 

critical  field 

+.14  ±  3.08 

4.46 

Harris  (1969) 

critical  field 

+.41  ±  4.72 

4.73 

0.68* 

*free  electron  gas  value. 


There  appears  to  be  a  great  deal  of  scatter  in  the  data  for  dlnN  /dp. 
The  value  of  Garfinkel  &  Mapother  is  apparently  not  correct  as  has 
been  established  by  private  communication  and  is  therefore  not  included 
in  Table  13.  The  data  obtained  from  critical  field  measurements  has 
a  number  of  limitations  largely  involved  with  the  extrapolations  of 
experimental  data  to  zero  temperature.  These  limitations  have  been 
alluded  to  by  Carbotte  (1971). 

To  the  present  time,  no  knowledge  of  the  pressure  dependence 
of  the  band  structure  density  of  states  exists  either  theoretical 
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or  experimental.  In  view  of  this,  the  usual  procedure  is  to  assume 
the  free  electron  value  of  dlnN^/dp  =  K^/3  with  no  real  justification. 

It  can  be  seen  from  column  four  of  Table  13  that  experimental 
values  for  dlnN^/dp  differ  considerably  from  those  of  the  free 
electron  case. 

Information  about  the  band  structure  density  of  states  can 
be  obtained  from  De  Haas  -  van  Alphen  measurements  and  such  data 
exists  for  lead  due  to  Anderson  &  Gold  (1965).  The  pressure  dependence 
of  the  Fermi  surface  has  been  studied  for  lead  by  Anderson  et  al 
(1966)  also  using  the  De  Haas  -  van  Alphen  effect.  They  have 
determined  the  pressure  dependence  of  the  Fermi  surface  in  lead  from 
which  it  is  possible  to  obtain  information  on  the  pressure  dependence 
of  the  band  structure  density  of  states  although  they  did  not  do  this. 
It  is  next  proposed  to  calculate  the  pressure  dependence  of  the  band 
structure  density  of  states  assuming  it  to  be  a  function  of  energy 
and  pressure  only,  thus:- 

Ng  =  NB(e,P)  6-26 


and  where  all  the  N's  used  are  per  unit  volume. 


From  6-26 

dlnNj 

D 

'91nN  ' 

D 

dEF  + 

'31nN  ' 

D 

l  dP  J 

F.S. 

l  de  J 

P  dP 

[  9P 

where  is  the  Fermi  energy. 

The  two  factors  in  the  first  terms  on  the  R.H.S.  of  equation 

6-27  are  known  experimentally  from  the  work  of  Anderson  and  Gold 

and  Anderson  et  al .  An  attempt  is  now  made  to  calculate  the  second 

term  (which  is  =  0  for  free  electrons). 

We  must  have  that:- 


(where  N/V  is  the  number  of  electrons  per  unit  volume  and  we  also 
know  that  the  number  of  electrons,  N,  must  be  conserved  upon 
decreasing  V,  i.e.  increasing  pressure).  The  experimental  evidence 
of  Anderson  &  Gold  however  shows  that:- 

<  0  6-29 
dp 

Thus  the  integral  on  the  left  hand  side  of  equation  6-28  would  decreas 
unless  Ng  changes  with  pressure  to  accommodate  the  situation.  This 
requires  to  increase  with  pressure.  Let  us  assume  that  the  density 
of  states  increases  then  according  to:— 

Vp)  =  a  V0)  6-30 


where  a  is  to  be  determined  and  the  numbers  in  parentheses  indicate 
the  pressure.  Thus:- 


and 


£F(P) 

aNB(0)de 

r  £F(0) 
NB(0)d£ 


=  _J_ 
V(P) 


6-31 


6-32 


since  V(P)  =  V  + 


and 


1  _3V  =  - 

V  3p 


we  have  that  V(P)  =  V[l-K^dpJ 

or  by  expanding  1  -  1_  [1  +  K  dp] 

V(P)  v  r 


equation  6-31  thus  becomes:- 
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•  Vp) 

aNB(0)de  =  N  +  N  I^dp 
V  V  1 


6-33 


Equation  6-32  may  be  written  in  the  form:- 


V0) 

*B(0)di 


The  last  term  in  equation  6-34  is  to  a  good  approximation 
,  eF(0) 


F 

f  F  r 

NB(0)d£  = 

NB(0)d£  + 

'  0  J 

e 

6-34 


NB(0)de  =  Np  dev  dp 


£F(P) 


dp 


where  Np  is  the  average  density  of  states  at  the  Fermi  surface, 
Thus  equation  6-34  becomes:- 


£F(0> 
NB(0)de  = 


.  eF(P) 

Ng (0)de  -  Np  d£p  dp  =  N 
dp  V 


6-35 


Substituting  this  expression  in  equation  6-33  for  N/V  gives:- 


eF(P) 

aN„(0)de  = 

D 


■  eF(p) 

NB(0)de  -  Np  deF  dp  +  N  Kpdp 


u 

Vp> 

up 

Ep(P)  ,1 

NB(0)d£  -  Np  d£p  dp  +  Kp 

dp 

NB(0)de  +  0 

o 

dp 

_  o  J 

where  we  have  used  equation  6-35  for  N/V.  Neglecting  the  second 
order  small  quantity  gives :- 


[a-l-Kpdp] 


eF(P) 

NB(0)de  =  -Np  deF  dp 


6-36 


dp 


making  the  reasonable  assumption  that 

i 

6-36  and  re-arranging  terms  gives :- 
a  =  1  +  dpfKT  -  Np  d£p 


eF(P) 

NB(0)de 


N  in  equation 
V 


N/V  dp 


6-37 
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Now  by  differentiating  equation  6-30  we  get:- 


dlnNB(P)  =  dot 


dp 


dp 


6-38 


which  using  equation  6-37  gives :- 

dlnNB(P)  =  K.J,  -  Np  dep  6-39 

~dp  N Tv  dp“ 

We  have  now  evaluated  a  form  for  the  second  term  in  equation  6-27 
so  that  the  final  expression  obtained  for  the  pressure  dependence 
of  the  band  structure  density  of  states  at  the  Fermi  surface  is:- 


dlnN 


B 


dp 

\  r  J 


F.S. 


dlnN 


B 


de 


-  N, 


N/V 


deF  + 


dp 


or 

dlnN  ' 

D 

= 

£F 

dlnN^ 

D 

-  e  N 

F  F 

dineF 

dp 

F.S. 

[  de  J 

P  N'v 

dp 

6-40 


The  validity  of  this  relationship  can  be  tested  for  the  case  of  the 
free  electron  gas.  In  this  case  the  terms  on  the  right  hand  side  of 
equation  6-40  are:- 


dlnN 


B 


de 


-  u 


6-41 


f*  =  4 

N/V  Z 


dln£p  =  2l  K^. 


6-42 


6-43 


dp 


and  upon  substitution  yield :- 


dlnN 


B 


=  .68  x  10  ^  bar 


dP  ) 


F.S. 


which  is  in  complete  agreement  with  the  value  of  Table  13. 

Turning  now  to  the  case  of  lead,  the  work  of  Anderson  &  Gold  shows 
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the  term  given  by  equation  6-41  to  be  — 2.5  times  that  for  the  free 
electron  case.  The  second  term  (given  by  equation  6-42)  is 
approximately  the  same  as  that  for  free  electrons  while  the  third 
term  (equation  6-43)  has  been  determined  by  Anderson  et  al  and 
found  to  be  (dln£^/<ip) 


f  '~t'/F.S. 


=  -.44  x  10  6  bar  1 


Using  these  values  in  equation  6-40  gives  for  pure  lead:— 


dlnN 


B 


dp 

\  L  / 


=  3.26  x  10  6  bar  1 


6-44 


F.S. 


It  can  be  seen  from  Table  13  that  this  value  is  in  very  good 
agreement  with  the  experimental  values  found  in  column  four. 

Using  the  definition  of  the  heat  capacity  coefficient  y  given 
by 

2  2 

Y  -  2tt  N^Zk^  (per  unit  volume) 


it  can  be  shown  that 

dlny  =  dlnN  +  dlnZ  at  the  Fermi  surface 
dP  dP 

By  comparing  this  with  equation  6-22  we  see  that 
dlnN  =  dlny 


6-45 


dp 


dp 


Hansen  et  al  (for  their  12%  data)  calculate  a  value  of  dlny  using 

dp 

equation  6-45  and  using  a  value  of  dlnN-  for  free  electrons  of  -1. 

D  — 

dlnV  3 

They  then  erroneously  compare  their  value  of  dlny  =  4.0  ±  1  x  10  bar 

dp 

with  that  of  White  and  claim  good  agreement.  The  comparison  however 
is  not  valid  since  White  assumes  the  volume  dependence  of  the  band 


137 


structure  density  of  states  to  be  +2/3  for  the  free  electron  case 
thus  indicating  the  explicit  inclusion  of  the  volume  in  his 
calculations  as  opposed  to  those  of  Hansen  et  al  whose  calculated 
value  is  per  unit  volume. 

Their  value  should  therefore  be  compared  with  the  value 
appearing  in  Table  13  in  which  case  the  agreement  is  poor. 

The  inescapable  conclusion  is  that  the  use  of  the  free  electron 
value  for  the  band  structure  density  of  states  is  not  really 
justified  and  that  it  should  be  replaced}  for  lead  at  least,  by  a 
value  such  as  derived  in  equation  6-44  which  supports  the 
experimental  observations. 


CHAPTER  VII 


CONCLUSIONS 


In  yiew  of  the  fact  that  conclusions  on  the.  various  results 
obtained  permeate  the  discussions  of  Chapter  VI,  only  the  salient 
features  are  repeated  here. 

The  PbxIn^_^_  system  exhibits  a  similar  trend  to  weaker  coupling 
as  a  function  of  applied  pressure  as  can  be  seen  from  several  view¬ 
points.  The  observed  reduction  in  the  gap  ratio  for  example 
illustrates  this  directly. 

The  phonon  spectra  and  gap  function  obtained  from  the  inversion 
process  indicate  a  general  increase  to  higher  energies  as  well  as  a 
decrease  in  structure  amplitude  (although  the  effect  is  less 
pronounced  than  for  the  samples  containing  lesser  concentrations  of 
indium).  This  is  in  good  agreement  with  the  theoretical  prediction 
of  Trofimenkoff  and  Carbotte  who,  using  a  scaled  model  phonon  spectrum 
similar  to  those  that  are  actually  observed,  predict  a  value  for 
dlnA/dlnTc  -  2.0  compared  with  the  values  of  1.38  ±.53  and  1.52  ±.58 
for  the  Pbggln^  and  Pb^^In^g  alloys  respectively.  The  strength 
electron-phonon  interaction  parameter  X  was  itself  observed  to 
decrease  indicative  of  the  general  trend  to  weaker  coupling. 

The  Gruneisen  parameters  calculated  are  in  very  good  agreement 
with  those  of  other  workers.  It  is  noticeable  that  with  increased 
indium  concentration  the  Gruneisen  y  associated  with  the  localised 
mode  decreased  considerably.  This  is  in  good  agreement  with  the 
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findings  of  Galkin  et  al  and  is  possibly  tied  in  with  the  findings 
of  Sood  which  show  the  local  mode  of  vibration  to  exhibit  at  least 
two  peaks  in  the  phonon  spectrum  and  whose  structure  and  location 
are  concentration  dependent. 

A  value  for  the  pressure  dependence  of  the  band  structure 
density  of  states  at  the  Fermi  surface  has  been  derived,  and  was 
found  to  be:- 

dlnNfi  =  3.26  x  10  ^  bar  ^  for  lead 
dp 

This  value  is  suggested  as  a  replacement  for  the  free  electron  value 
“•6  —I 

of  0.68  x  10  bar  in  future  theoretical  considerations  of  lead  and 
lead  like  systems. 

Throughout  the  experiments  performed  the  tunnel  junctions  have 
themselves  provided  the  most  reliable  method  of  pressure  determination 
both  in  the  solid  and  in  the  fluid.  Their  use  in  the  solid  helium 
bath  provided  direct  evidence  that  this  technique  of  pressure 
application  does  indeed  result  in  hydrostatic  conditions. 
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